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Abstract — The design of modulation schemes for the physi- 
cal layer network-coded two way wireless relaying scenario is 
considered. It was observed by Koike-Akino et al. |4| for the 
two way relaying scenario, that adaptively changing the network 
coding map used at the relay according to the channel conditions 
greatly reduces the impact of multiple access interference which 
occurs at the relay during the MA Phase and all these network 
coding maps should satisfy a requirement called exclusive law. We 
extend this approach to an Accumulate-Compute and Forward 
protocol which employs two phases: Multiple Access (MA) phase 
consisting of two channel uses with independent messages in each 
channel use, and Broadcast (BC) phase having one channel use. 
Assuming that the two users transmit points from the same 4- 
PSK constellation, every such network coding map that satisfies 
the exclusive law can be represented by a Latin Square with 
side 16, and conversely, this relationship can be used to get the 
network coding maps satisfying the exclusive law. Two methods 
of obtaining this network coding map to be used at the relay are 
discussed. Using the structural properties of the Latin Squares 
for a given set of parameters, the problem of finding all the 
required maps is reduced to finding a small set of maps. Having 
obtained all the Latin Squares, the set of all possible channel 
realizations is quantized, depending on which one of the Latin 
Squares obtained optimizes the performance. The quantization 
thus obtained, is shown to be the same as the one obtained in 
|7| for the 2-stage bidirectional relaying. 

I. Background 

The concept of physical layer network coding has 
attracted a lot of attention in recent times. The idea of 
physical layer network coding for the two way relay channel 
was first introduced in HI, where the multiple access 
interference occurring at the relay was exploited so that the 
communication between the end nodes can be done using 
a two stage protocol. Information theoretic studies for the 
physical layer network coding scenario were reported in |2), 
0. The design principles governing the choice of modulation 
schemes to be used at the nodes for uncoded transmission 
were studied in J4j. An extension for the case when the 
nodes use convolutional codes was done in [5]. A multi-level 
coding scheme for the two-way relaying was proposed in [6]. 

We consider the two-way wireless relaying scenario shown 
in Fig. 1, where two-way data transfer takes place among the 
nodes A and B with the help of the relay R. It is assumed 
that the two nodes operate in half-duplex mode, i.e., they 
cannot transmit and receive at the same time in the same 
frequency band. The relaying protocol consists of two phases, 



Fig. 1. A two-way ACF relay channel 

multiple access (MA) phase, consisting of two channel uses 
during which A and B transmit to R twice, two independent 
messages in the two channel uses, with points from 4-PSK 
constellation, and broadcast (BC) phase, in which R transmits 
to A and B. The relay node R accumulates the information 
sent by the user nodes in the first and second channel use of 
the MA phase, and transmits in the BC phase a message that 
contains information about all the four messages received by 
it in the MA phase. Network Coding is employed at R in such 
a way that A(/B) can decode the two messages ttansmitted 
by B(/A), given that A(/B) knows its own messages. We call 
this strategy accumulate-compute and forward (ACF) protocol. 

It was observed in |4| and [7| for 4-PSK, that for uncoded 
transmission, the network coding map used at the relay 
needs to be changed adaptively according to the channel fade 
coefficient, in order to minimize the impact of multiple access 
interference. In other words, the set of all possible channel 
realizations is quantized into a finite number of regions, with 
a specific network coding map giving the best performance 
in a particular region. It is shown in JS] for any choice of 
signal sets of equal cardinality used at the two users, that 
every such network coding map that satisfies the exclusive 
law is representable as a Latin Square and conversely, this 
relationship can be used to get the network coding maps 
satisfying the exclusive law. 

Definition 1: A Latin Square of order M is an M x M 
array in which each cell contains a symbol from a set of t 
different symbols such that each symbol occurs at most once 
in each row and column (9). 

Similar to the ACF protocol, a store-and-forward protocol 
has been earlier studied in [10|, for the two-way relaying 
channel. In IflOl , the authors derive an upper bound on the 
ergodic sum-capacity for the two-way relaying scenario when 
delay tends to infinity, and propose two alternative awaiting 
and broadcast (AAB) schemes which approach the new upper 
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F 2 XF 2 = -| 00, 01, 10, 11} 
Z4= {0, 1, 2, 3 } 

Fig. 2. 4-PSK constellation 

bound at high SNR. Using numerical results, they show that 
the proposed AAB schemes significantly outperforms the 
traditional physical layer network coding methods without 
delay in terms of ergodic maximum sum rates. However, 
modulation and physical layer network coding have not been 
addressed in iflOl . 

The remaining content is organized as follows: Section II 
discusses the basic concepts, definitions and a summary of 
the contributions of this paper. Section III demonstrates the 
network code obtained using Cartesian Product that is utilized 
at the relay for two-way ACF relaying which removes the 
fade states associated with the channels. In Section IV, we 
show how this network code can be obtained using Singularity 
Removal Constraints. Section V gives results based on struc- 
tural properties of Latin Squares. In Section VI the complex 
plane is quantized depending on which one of the obtained 
Latin Squares maximizes the minimum cluster distance and 
Section VII gives the simulation results that demonstrate the 
improvement in the performance using the suggested scheme. 
Section VIII concludes the paper. 

II. Preliminaries 

Let S denote the symmetric 4-PSK constellation {±1 ± 1} 
as shown in Fig. 2, used at A and B. Assume that 
A(/B) wants to send two 2-bit binary tuples to B(/A). 
Let \i : F 2 , — > S denote the mapping from bits to 
complex symbols used at A and B where F2 = {0, 1}. Let 
xa x = fi(sA ± ) ,XBi = ^(ssj € S denote the complex 
symbols transmitted by A and B at the first channel use 
respectively, and xa 2 — n( s A 2 ) i x b 2 = ^{ s b 2 ) £ S denote 
the complex symbols transmitted by A and B at the second 
channel use respectively, where sa 11 sb 11 sa 21 sb 2 & F 2 ,. 

Multiple Access (MA) Phase: 

It is assumed that the channel state information is not 
available at the transmitting nodes A and B during the MA 
phase. The received signal at R at first channel use is given 
by 

Y Rl = H A x Al + H B x Bl + Z Rl (1) 
and the received signal at R at the second channel use, 

Yr 2 = H a xa 2 + H b xb 2 + Z R2 (2) 

where Ha and Hr are the fading coefficients associated with 
the A-R and B-R link respectively. Note that we are taking 



Ha and Hb to be the same for the two channel uses. The 
additive noise Z Rl and Zr 2 are assumed to be CAT (0,<7 2 ), 
where CM (0, cr 2 ) denotes the circularly symmetric complex 
Gaussian random variable with variance a 2 . We assume a 
block fading scenario, with z — ^e^ 9 = Hb/Ha, where 
7 G R + and — n < 9 < ir, is referred to as the fade state for 
the first and second transmission by A and B at the first and 
second channel use, and for simplicity can also be denoted 
by (7, 9). Also, it is assumed that z is distributed according 
to a continuous probability distribution. 

Let Sr (7, 9) denote the effective constellation seen at the 
relay during the MA phase, i.e., 

Sr(pi,6) = {(xi +-ye 39 y i ,x j +je je y j )\x i) yi,Xj,yj e S} . 

The effective constellation remains the same over the two 
channel uses, since we assume Ha and Hb and hence the 
ratio Hb /Ha = ^e? to be the same during the two channel 
uses. 

Let d m in (j^ 9 ) denote the minimum distance between the 
points in the constellation Sr (7, 9) during MA phase, as given 
by © on the next page. From (01, it is clear that there exist 
values of je^ 9 , for which d m i n (je^ 9 ) = 0. Let, 

H = {jeJ 9 e C\d mm ( 7 e*) - 0} . 

The elements of H are called singular fade states. For singular 
fade states, |«Sr ( 7 ,0)| < 4 4 . 

Definition 2: A fade state ^e^ 9 is defined to be a singular 
fade state for the ACF two-way relaying, if the cardinality of 
the signal set Sr (7, 9) is less than 4 4 . Let % denote the set 
of singular fade states for the two-way ACF relaying. 

Let (xAi , %Bx ) and (x a 2 , xb 2 ) € S 2 denote the Maximum 
Likelihood (ML) estimate of (lAn^Bi) and (xa 2 ,xb 2 ) at 
R based on the received complex numbers Yr x and Yr 2 at 
the two channel uses, as given in 

Broadcast (BC) Phase: 

Depending on the value of je^ 6 , R chooses a map 
M. 1 ' 9 : 6> 4 — > S' where S is a complex signal set of size 
between 4 2 and 4 4 used by R during the BC phase. 

The received signals at A and B during the BC phase are 
respectively given by, 

Y A = H' A X R + Z A and Y B = H' B X R + Z B (10) 

where Xr = M 1 ' 9 ((i^ , xbJ , (xa 2 ,xb 2 )) € S' is the 
complex number transmitted by R. The fading coefficients 
corresponding to the R-A and R-B links are given by H' A 
and H' B respectively and the additive noises Za and Zb are 
CM(0,a 2 ). 

The elements in 6> 4 which are mapped to the same signal 
point in S' by the map M. 1 ' 9 are said to form a cluster. 
Let {£1, £2, ■•, £1} denote the set of all such clusters. The 
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d 2 m l n(^ jS )= min j I (x Al -x' A ) + 'yeP (x Bl ~ z' s ) 1 2 + I (xa 2 - ) + -ye> e (. 

(("A! .a:B 1 )>(^A2 ' a: - B 2 D^tt^! ^Bj )<( X A 3 > x B 2 )) 
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((iipiBj), (xa 2 ,x b )) = arg min { W R - H A x' A - H B x' B I + |y H - Ha^a' - h bx' b ,\ \ (4) 
((^ 1 ^B 1 M^ 2 .^B 2 ))e.S 4 I 1 11 2 1 J 

.M 7 ' 9 ((xa 1 ,x j 4 2 ) , (iBpiBj)) A4 7 ' 9 ((^Aj.^As) ' ( x Si,a:B 2 )) , uj/ieneuer (ia!,^ 2 ) ^ ( I l4 1 i I A 2 ) Vx Bi ,xb 2 G 5 (5) 

■M 7 ' 9 ((xa,,^) , (^ Bi ,^b 2 )) ^ -M 7,9 ((iAi^A 2 ) i (^ Bi ^b 2 )) i whenever (x Bi ,xb 2 ) ¥= {^'b 1 ' x 'b 2 ) ^x Al ,x A . 2 e S (6) 

ss? (v)) 2 = ^ zn <s - ^i) + 7e - r + 1 k - ^ - *oTf z~ 

V V // ((lA 1 ,I Bl ))(^ ! ^fl 2 ))6£iU V 7 V 71 IV 7 V 71 J 
<X*Ai • x 'b 1 )'( x 'a 2 ' x 'b 2 )) < ^ c j 

dLn (C 1 ' 6 ) = mki " " " { I (x Al - S + 7e jS (x Bl - s' fll ) I* + I (xa 2 - x' A ) + ^ (xb 2 - ^ ) P I 



((SA] .»B 1 ),( a: A2 • X B 2 ))-(( X A 1 • x B 1 )>( x 'a 2 > x 'b 2 ))£ S y 
M^- e ((x Al ,x Bl ),(xA 2 ,XB 2 ))^M~<- t> ax Al ,x' Bl ),(x A2 , X ' B2 )) 



(8) 



dLn(c h ,ie 3t> ) = min J (a;^ - x' A ) + -ye' 6 (x Bl - a>' B )\ + I (xa 2 - x'a 2 ) + ie i$ (x B2 - x' B 
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formation of clusters is called clustering, denoted by C ieJ . 

In order to ensure that A(/B) is able to decode B's(/A's) 
messages, the clustering C 1&1 should satisfy the exclusive 
law, as given in (|5), (JSJ above. 

Definition 3: The cluster distance between a pair of clusters 
Ci and Cj is the minimum among all the distances calculated 
between the points ((xa 1 ,Xb 1 ) , (xa 2 ,Xb 2 )) € £, and 
( (x' Ai , x' Bl ) , (x' A2 , x' B2 ) ) G £j in the effective constellation 
used by the relay node R, as given in (0 above. 

Definition 4: The minimum cluster distance of the 
clustering C 76 is the minimum among all the cluster 
distances, as given in dHJ above. 

The minimum cluster distance determines the performance 
during the MA phase of relaying. The performance during 
the BC phase is determined by the minimum distance of the 
signal set S . For values of je^ 6 in the neighborhood of the 
singular fade states, the value of d m i n (c ie3 ^ is greatly 
reduced, a phenomenon referred to as distance shortening [4]. 
To avoid distance shortening, for each singular fade state, a 
clustering needs to be chosen such that the minimum cluster 
distance is non zero. 

A clustering C h is said to remove singular fade state 
h € if d m in (C ) > 0. For a singular fade state h € %, 
let C h denote the clustering which removes the singular fade 
state h (if there are multiple clusterings which remove the 
same singular fade state h, choose any of the clusterings). 
Let C n = {C h : h € Ti} denote the set of all such clusterings. 

Definition 5: The minimum cluster distance of the cluster- 
ing C h ,h G H at the fade state -fe^ 9 which is not necessarily 
a singular fade state, denoted by d m i n {C h , r ye^ s '), is as given 



in ©. Note that if -ye 10 = h G H, d min (C h ,h) , reduces to 

u >nin 

(C h ) given in ©. 

In general, the channel fade state ^d* need not be a singular 
fade state. In such a scenario, among all the clusterings 
which remove the singular fade states, the one which has the 
maximum value of the minimum cluster distance at 7e- J is 
chosen by the relay R. In other words, for 7e j9 ^ W, the 
clustering C 7 ' 61 is chosen to be C h by the relay R, which 
satisfies d mm (C h ,^ e ) > d mm (c h ' ,je^ ,Vft + hi e H. 
Since the clusterings which remove the singular fade states 
are known to all the three nodes and are finite in number, the 
clustering used for a particular realization of the fade state can 
be indicated by R to A and B using overhead bits. 

In (8), such clusterings that remove singular fade states 
for two-way 2-stage relaying scenario were obtained with the 
help of Latin Squares. 

The contributions of this paper are as follows: 

• It is shown that if the users A and B transmit points from 
the same 4-PSK constellation, the clusterings proposed in 
|8| for two user case can be utilized to get clusterings for 
this case for removing the singular fade states containing 
either 16 or 25 points by introducing the notion of 
Cartesian Product of Clusters. In other words, the 16 x 16 
Latin Squares representing the ACF relaying can be 
obtained with the help of 4 x 4 Latin Squares representing 
the clusterings for two-way 2-stage relaying as given in 

m. 

• Another clustering is proposed for the ACF protocol 
in the two-way relay channel called Direct Clustering. 
This clustering also removes the singular fade states and 
reduces the number of clusters for some cases. Using this 
clustering, the size of the resulting constellation used by 
the relay node R in the BC phase is reduced to 20 for a 
category of cases, as compared to the Cartesian Product 
approach which results in the constellation size being 25 
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for these cases. 

• The quantization of the complex plane that contains all 
the possible fade states, depending on which one of 
the obtained clusterings maximizes the minimum cluster 
distance, is proven to be the same as for the two-way 
2-stage relaying scenario as done in Q. 

• Simulation results indicate that at high SNR, the schemes 
based on the ACF protocol performs better than the 
schemes proposed in (4|, based on two-stage two 
way relaying. With 4-PSK signal set used at the end 
nodes, the ACF protocol achieves a maximum sum 
throughput of 8/3 bits/s/Hz, whereas it is 2 bits/s/Hz for 
the schemes based on 2-stage two way relaying. 



III. Exclusive Law and Latin Squares 

The nodes A and B transmit symbols from the same 
constellation, viz., 4-PSK. Our aim is to find the map that the 
relay node R should use in order to cluster the 4 4 possibilities 
of ((iaj, xb ± ) , (xa 2 , x b 2 )) sucn that the exclusive law given 
by ©, © is satisfied. Consider the 16 x 16 array consisting 
of the 16 possibilities of (xa 1 ,%a 2 ) along the rows and the 
16 possibilities of (xb 1 ,xb 2 ) along the columns. We fill this 
array with elements from C — {C±, £2, ■ £t} where each 
symbol denotes a unique cluster. The constellation size used 
by the relay in the BC phase, or the number of clusters has 
to be at least 16, since each user needs 4 bit information 
corresponding to the two messages sent by the other user 
implying t > 16. In order to keep (O and © satisfied, a 
symbol from C can occur at most once in each row and 
column. So the 16 x 16 array having (xa i: xa 2 ) along the 
rows and (xb x ,xb 2 ) along the columns must be a Latin 
Square of side 16 (Definition 1). The equivalence between 
the network code used by the relay in two-way relaying 
scenario and Latin Squares has been previously discussed in 
[7]. The clusters are obtained by putting together all those 
((xAt , xb 1 ) , (xa 2 , xb 2 )) for which the corresponding entry 
((xAt, x a 2 ), (xb 1 ,Xb 2 )) in the array is the same. 

From above, we can say that all the relay clusterings that 
satisfy the mutually exclusive law forms Latin Squares of 
order 16 with entries from C with t > 16, when the end 
nodes use PSK constellations of size 4. It therefore suffices 
to consider the network code used by the relay node in the 
BC phase to be a 16 x 16 array with rows(/columns) indexed 
by the 2-tuple consisting of the symbols sent by A(/B) during 
the first and second channel use. The cells of the array must 
be filled with elements of C in such a way, that the resulting 
array is a Latin Square of order 16 and t > 16. 

Removing Singular fade states and Constrained Latin 
Squares 

The relay can manage with constellations of size 16 in 
BC phase, but it is observed that in some cases relay may 
not be able to remove the singular fade states and results in 
severe performance degradation in the MA phase. As stated 
in Section II, that a clustering C h is said to remove singular 



fade state h € %, if d m i n (C h ) > 0. Removing singular fade 
states for a two-way ACF relay channel can also be defined 
as follows: 

Definition 6: A clustering C h is said to remove 
the singular fade state h € %, if any two 
possibilities of the messages sent by the users 

((x Al ,x Bl ) ,(xa 2 ,xb 2 )) ,((x' Al ,x' Bi ) ,(x' A2 ,x' B2 )) e S 4 
that satisfy 

, _ x' Ai — X Al _ x' A2 — XA 2 
XB\ x b 1 Xb 2 Xg 2 

are placed together in the same cluster by C h . 

Definition 7: A set {((x Al , Xb x ) , (xa 2 , x b 2 ))} Q S 4 
consisting of all the possibilities of ((xa 1 ,xb 1 ),(xa 2 ,xb 2 )) 
that must be placed in the same cluster of the clustering 
used at relay node R in the BC phase in order to remove the 
singular fade state h is referred to as a Singularity Removal 
Constraint for the singular fade state h in two-way ACF 
relaying scenario. 

As given in (4), a complex number ^e^ 6 is 
defined to be a singular fade state for the two- 
way 2-stage relaying scenario, if xa + 7e J Xb = 
x' A + rye^x'g for some (xa, xb), (x' a , x' b ) £ S 2 and 
the set {(xa,xb) \ (xa,xb)&S 2 } consisting of all the 
possibilities of (xa,xb) G S 2 that must be placed in the same 
cluster of the clustering that removes the fade state je^ 6 for 
the two-way 2-stage relaying scenario is the corresponding 
set of singularity removal constraint for the singular fade 
state je^ 8 . As we show in the following lemma, the singular 
fade state for the ACF two-way relaying are the same as the 
singular fade state for the two-way 2-stage relaying scenario. 

Lemma 1: The singular fade states for the ACF two-way 
relaying scenario are the same as the 12 singular fade states 
for two-way 2-stage relaying scenario as computed in |0J. 

Proof: Let 'ye^ 6 be a singular fade state for 
the ACF two-way relaying scenario. By definition, 

3((xa 1 ,xb 1 ),(xa 2 ,xb 2 )),((x'a 1 ,x' Bi ),(x'a 2 ,x i b 2 )) e S 4 
such that, 

XA! +je je x Bl =x' Ai +je je x Bi , and 
xa 2 + ie je xB 2 = x' M + ~fe j9 x' B2 . 

where x Al , xb x , xa 2 , xb 2 , x' Al , x' Bi , x\ 2 , x' B<2 G S. 

Then, by definition, ^e^ 6 must be a singular fade state for 
two-user 2-stage relaying. 

Conversely, let ^e^ 6 be a singular fade state for two-user 
2-stage relaying scenario. Then, 3(2,4, xb), (x' a ,x' b ) £ S 2 
such that, 

x A + ^eJ e x B = x' A + 7c 7 x' B . 

Then, since for any (x, y) € S, x + je^y = 
x + je j6 y, {((xA,x B ),(x,y)), ((x' A ,x' B ), (x,y))} is a 
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subset of a singularity removal constraint and je^ 9 is a 
singular fade state for the ACF two-way relaying scenario. 

Thus, the singular fade states for the ACF two-way relaying 
scenario and the singular fade states for the two-way 2-stage 
relaying scenario are the same. ■ 

Let je^ 9 be a fade state for the two-way ACF relaying 
scenario. Then je^ 9 can be viewed as a fade state for the 
first and second channel use in the MA phase as shown in 
Lemma 1. In 0) and Q, it is shown that for the two-way 
2-stage relaying, the 4 2 possible pairs of symbols from 
4-PSK constellation sent by the two users in the MA phase, 
can be clustered into a clustering dependent on a singular 
fade coefficient, of size 4 or 5 in a manner so as to remove 
this singular fade coefficient. In the case of two-way ACF 
relaying, at the end of MA phase, relay receives two complex 
numbers, given by ([TJ and Instead of R transmitting a 
point from the 4 4 point constellation resulting from all the 
possibilities of ((x^, xb x ) , (xa 2 i x b 2 )) me relay R can 
choose to group these possibilities into clusters represented by 
a smaller constellation. One such clustering for the case when 
7e J can be obtained by utilizing the clustering provided in 
[7] for the two-way 2-stage problem in order to remove this 
fade state. Let denote the clustering for the physical 
network coded two-way relaying scenario that removes the 
singular fade state h E C for the two-way 2-stage relaying 
case as given in [7]. 

Definition 8: We define the Cartesian Product of a cluster- 
ing = {h,l 2 , -,l m } with itself denoted by V^ h \ where 
for i = 1,2, m; 

h { {Xi 1 , Ui 1 ) , (Xi 2 , Ui 2 ) , (%i Bi j Vi Bi ) } 

with x ip ,y ip eZ 4 Vp = 1,2, ...,Sj as follows: 

v [h] _ ^c {luh} , ...,C {h ' lm} , ...,C {lm ' h} , ■■■,C {lm ' lm} } 
where, 

ci l *M = {{(x ip ,y ip ),(x jq ,y jg )) I p = 1,2,.., Si and q = 1, 2, .., Sj } . 

Lemma 2: Let 7e je € H. The clustering obtained by taking 
the Cartesian Product T>b ei$ ] of with itself removes 

the singular fade state ^e' for the two-way ACF relaying 
scenario. 

Proof: Let C^] = Z 2 , ... ,Z m },where for i = 
1,2,..., m, 

h { {Xi 1 , Ui 1 ) , (x*i 2 , Ui 2 ) , (Xi B , Vi Si ) } 

Then, 

vb eSe ] = |c {h -' i} , ...,c {h ' lm} ,...,c {lmM} , ...,c {lm ' lm A 

where, 

Cih'lj} = {(( Xip ,y ip ),(x jq ,y jq )) \p = l,2,.., Si and q = 1, 2, .., Sj } . 

By definition, a singularity removal constraint for the fade 
state 7e J in the case of two-way ACF relaying scenario is a 



set {((xi,yi), | i = 1, 2, *} such that VI < i 1: i 2 < 

t, 

id x ii ~ x ii X ii ~~ X ii 

l e = ~ = / _ / 

Vi\ Ui2 Vi 1 Vi 2 

Now, 7e j9 = l^Z^ 1 => (xi,yi) and (x 2 ,2/ 2 ) must belong 
to the same cluster, say k in cb el J for it to remove the fade 

state 7e j6 '. Similarly, -fe^ 9 = rE ' 2 _ z ' 1 =>■ (cc' x , ) and (a; 2 ,y 2 ) 

y%i y% 2 

must belong to the same cluster, say lj in d ie3 "]. This holds 
VI < i% < t. Thus, the singularity removal constraint is, 

{((n.IttMzi.i/i)) M = i,2,...,t}cc^^>, 

for some 1 < i,j < n. Therefore, the clustering T>b e3 ] 
removes the singular fade state -fe^ 9 . ■ 

It was studied in [4] by computer search, and then in [7] 
analytically, that there are 12 possible singular fade states 
in the complex plane for the case when two users transmit 
points from 4-PSK constellation in the MA phase. Out of 
these 12 fade states, 4 lie on the unit circle, 4 lie on a circle 
of radius \J~2 and 4 lie on a circle of radius l/y2. The size of 
the constellation used at R in the BC phase for these cases is 
either 4 or 5. For the two-user ACF scenario we are dealing 
with, we have two channel uses in the MA phase, with A 
and B transmitting a message each in the first channel and 
second channel uses, with the two messages sent by each user 
in the first and second channel use being possibly different. 
Keeping in mind the three classes of fade states depending on 
the radius of the circle it lies on, we consider the following 
three cases: 

Case 1: je^ 9 lies on the unit circle. 

Case 2: ^e> 9 lies on the circle of radius l/y/2. 

Case 3: ^e> 9 lies on the circle of radius \/2. 

Case 1: je j9 lies on the unit circle. 

Since both user nodes A and B require 4 bits of information 
from the other user, the size of the constellation that R uses 
will be at least 2 4 = 16. The Cartesian Product of C^ 3 "] 
with itself consists of 16 clusters since the clustering cb eJ ] 
has 4 clusters. We illustrate this case with the help of the 
following example, the remaining instances of fade states that 
lie on the unit circle can be obtained using this example, as 
we show in Lemma 6 in Section V. 

Example 1: Let the fade state 7e j6 = j. The clustering 
for the case as given in [ 8 1 is given by, 

C b] ={h,h,h,U} where, 

h = {(0,0), (1,2), (2,1), (3, 3)} 

1 2 = {(0,3), (1,1), (2, 2), (3,0)} 

1 3 = {(0,1), (1,3), (2,0), (3, 2)} 
k = {(0,2), (1,0), (2, 3), (3,1)}. 
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The Cartesian Product of the above clustering given by lM = 
ij} | ^ j = 1,2,3,4} contains exactly 16 clusters: 



{((0,0) 


(0,0)), 


((0,0), (1,2 


), 


((0,0), (2,1)), 


((0,0), (3,3) 


((1,2), 


(0,0)), ( 


(1,2) 


(1,2) 


), ( 


(1,2), 


(2, 1)), ( 


(1, 2) 


(3,3)) 


((2, 1), 


(0,0)), ( 


(2, 1) 


(1, 2) 


• \ 


(2,1), 


(2, 1)), ( 


(2, 1) 


(3,3)) 


((3,3), 


(0,0)), ( 


(3, 3) 


(1, 2) 


),( 


(3, 3), 


2, 1)), ( 


(3, 3) 


(3, 3)) 


{((0,0) 


(0,3)), 


((0,0) 


,(1,1 


), 


((0,0) 


(2,2)), 


((0,0) 


,(3,0) 


((1,2), 


(0,3)), ( 


(1,2) 


(1, 1) 


• I 


(1,2), 


(2,2)), ( 


(1, 2) 


(3,0)) 


((2, 1), 


(0,3)), ( 


(2, 1) 


(1, 1) 




(2,1), 


(2,2)), ( 


(2, 1) 


(3,0)) 


((3,3), 


(0,3)), ( 


(3, 3) 


(1, 1) 


),( 


(3, 3), 


2, 2)), ( 


(3, 3) 


(3, 0)) 


{((0,0) 


(0,1)), 


((0,0) 


,(1,3 


), 


((0,0) 


(2,0)), 


((0,0) 


, (3, 2) 


((1, 2), 


(0, 1)), ( 


(1,2) 


(!,3) 




(1,2), 


(2,0)), ( 


(1, 2) 


(3,2)) 


((2, 1), 


(0, 1)), ( 


(2, 1) 


(!, 3) 




(2,1), 


(2,0)), ( 


(2, 1) 


(3, 2)) 


((3,3), 


(0, 1)), ( 


(3, 3) 


(1, 3) 






(2, 0)), ( 


(3, 3) 


(3, 2)) 


{((0,0) 


(0,2)), 


((0,0) 


,(1,0 




((0,0) 


(2,3)), 


((0,0) 


, (3, 1) 


((1,2), 


(0,2)), ( 


(1,2) 


(!,0) 




(1,2), 


(2,3)), ( 


(1, 2) 


(3, 1)) 


((2, 1), 


(0,2)), ( 


(2, 1) 


(!,0) 




(2,1), 


(2,3)), ( 


(2, 1) 


(3, 1)) 


((3,3), 


(0,2)),( 


(3, 3) 


(1, 0) 




(3,3), 


2, 3)), ( 


(3, 3) 


(3, 1)) 


{((0,3) 


(0,0)), 


((0,3) 


,(1,2 




((0,3) 


(2,1)), 


((0,3) 


,(3,3) 


((1, 1), 


(0,0)), ( 


(1, 1) 


(1,2) 




(1,1), 


(2, 1)), ( 


(1, 1) 


(3,3)) 


((2,2), 


(0,0)), ( 


(2, 2) 


(1, 2) 




(2,2), 


(2, 1)), ( 


(2, 2) 


(3,3)) 


((3,0), 


(0,0)), ( 


(3, 0) 


(1, 2) 




(3, 0), 


2, 1)), ( 


(3, 0) 


(3, 3)) 


{((0,3) 


(0,3)), 


((0,3) 


,(1,1 




((0,3) 


(2,2)), 


((0,3) 


,(3,0) 


((1, 1), 


(0,3)), ( 


(1, 1) 


(1, 1) 




(1,1), 


(2,2)), ( 


(1, 1) 


(3,0)) 


((2, 2), 


(0,3)), ( 


(2, 2) 


(1, 1) 




(2,2), 


(2, 2)), ( 


(2, 2) 


(3,0)) 


((3,0), 


(0,3)), ( 


(3, 0) 


(1, 1) 




(3,0), 


2, 2)), ( 


(3, 0) 


(3, 0)) 


{((0,3) 


(0,1)), 


((0,3) 


,(1,3 




((0,3) 


(2,0)), 


((0,3) 


,(3,2) 


((1, 1), 


(0, 1)), ( 


(1, 1) 


(!,3) 




(1,1), 


(2,0)), ( 


(1, 1) 


(3,2)) 


((2, 2), 


(0, !)),( 


(2, 2) 


(!,3) 




(2,2), 


(2,0)), ( 


(2, 2) 


(3,2)) 


((3,0), 


(0, !)),( 


(3, 0) 


(1, 3) 




(3, 0), 


(2, 0)), ( 


(3, 0) 


(3, 2)) 


{((0,3) 


(0,2)), 


((0,3) 


,(1,0 




((0,3) 


(2,3)), 


((0,3) 


, (3, 1) 


((1, 1), 


(0,2)), ( 


(1, 1) 


(!,0) 




(1,1), 


(2,3)), ( 


(1, 1) 


(3, 1)) 


((2,2), 


(0,2)), ( 


(2, 2) 


(!,0) 




(2,2), 


(2,3)), ( 


(2, 2) 


(3, 1)) 


((3,0), 


(0,2)),( 


(3, 0) 


(1, 0) 




(3,0), 


2, 3)), ( 


(3, 0) 


(3, 1)) 


{((0,1) 


(0,0)), 


((0,1) 


,(1,2 




((0, 1) 


(2,1)), 


((0,1) 


, (3, 3) 


((M), 


(0,0)), ( 


(!.3) 


(1,2) 




(1,3), 


(2, 1)), ( 


(1,3) 


(3,3)) 


((2,0), 


(0,0)), ( 


(2,0) 


(1, 2) 




(2,0), 


(2, 1)), ( 


(2,0) 


(3,3)) 


((3,2), 


(0,0)), ( 


(3, 2) 


(1, 2) 




(3,2), 


2, 1)), ( 


(3, 2) 


(3, 3)) 


{((0,1) 


(0,3)), 


((0,1) 


,(1,1 




((0, 1) 


(2,2)), 


((0,1) 


,(3,0) 


((M), 


(0,3)), ( 


(M) 


(1, 1) 




(1,3), 


(2,2)), ( 


(M) 


(3,0)) 


((2,0), 


(0,3)), ( 


(2,0) 


(1, 1) 




(2,0), 


(2,2)), ( 


(2,0) 


(3,0)) 


((3, 2), 


(0,3)), ( 


(3, 2) 


(1, 1) 




(3,2), 


2, 2)), ( 


(3, 2) 


(3, 0)) 


{((0,1) 


(0,1)), 


((0,1) 


,(1,3 




((0, 1) 


(2,0)), 


((0,1) 


, (3, 2) 


((M), 


(0, 1)), ( 


(M) 


(!,3) 




(1,3), 


(2,0)), ( 


(1,3) 


(3,2)) 


((2,0), 


(0, !)),( 


(2,0) 


(!,3) 




(2,0), 


(2,0)), ( 


(2,0) 


(3,2)) 


((3, 2), 


(0, 1)), ( 


(3,2) 


(!,3) 




(3,2), 


(2,0)), ( 


(3, 2) 


(3,2)) 


{((0,1) 


(0,2)), 


((0,1) 


,(1,0 




((0, 1) 


(2,3)), 


((0,1) 


, (3, 1) 


((M), 


(0,2)), ( 


(M) 


(!,0) 




(1,3), 


(2,3)), ( 


(1,3) 


(3, 1)) 


((2,0), 


(0,2)),( 


(2,0) 


(!,0) 




(2,0), 


(2,3)), ( 


(2,0) 


(3, 1)) 


((3, 2), 


(0,2)),( 


(3,2) 


(!,0) 




(3,2), 


(2,3)), ( 


(3,2) 


(3, 1)) 


{((0,2) 


(0,0)), 


((0,2) 


,(1,2 




((0, 2) 


(2,1)), 


((0,2) 


,(3,3) 


((!,0), 


(0,0)), ( 


(!,0) 


(1,2) 




(1,0), 


(2, 1)), ( 


(!,0) 


(3,3)) 


((2,3), 


(0,0)), ( 


(2,3) 


(1,2) 




(2,3), 


(2, 1)), ( 


(2,3) 


(3,3)) 


((3, 1), 


(0,0)), ( 


(3, 1) 


(1,2) 




(3,1), 


(2, 1)), ( 


(3, 1) 


(3,3)) 



= {((0,2), (0,3)) 
((1,0), (0,3)), 
((2, 3), (0,3)), 
((3, 1),(0,3)), 

= {((0,2), (0,1)) 
((1,0), (0,1)), 
((2, 3), (0,1)), 
((3, 1),(0,1)), 

= {((0,2), (0,2)) 
((1,0), (0,2)), 
((2, 3), (0,2)), 
((3, 1),(0,2)), 



,((0,2),(1 
((1,0), (1, 
((2,3),(1, 
((3,1), (1, 
,((0,2),(1 
((1,0), (1, 
((2,3),(1, 
((3,1),(1, 
,((0,2),(1 
((1,0), (1, 
((2,3),(1, 
((3,1),(1, 



>1)),((0,2) 
1)),((1,0), 
1)),((2,3), 
1)),((3,1), 
,3)), ((0,2) 
3)), ((1,0), 
3)), ((2,3), 
3)), ((3,1), 
0)),((0,2) 
0)),((1,0), 
0)),((2,3), 
0)),((3,1), 



• (2,2)),((0 

(2.2) ),((1 
(2, 2)), ((2 
(2, 2)), ((3, 
,(2,0)),((0 
(2,0)),((1 
(2,0)), ((2 
(2,0)), ((3, 
,(2,3)),((0 

(2.3) ),((1, 
(2, 3)), ((2, 
(2, 3)), ((3, 



2), (3,0)), 

0) ,(3,0)), 
3), (3,0)), 

1) ,(3,0))} 
2), (3, 2)), 

0) ,(3, 2)), 
3), (3, 2)), 

1) ,(3, 2))} 
2), (3,1)), 

0) ,(3, 1)), 
3), (3, 1)), 

1) ,(3, 1))} 



The entries of the above clusters are of the form 
((xa 1 ,xb 1 ),(xa 2 i x b 2 ))> i- e -> m the order A's transmission 
during the first channel use, B's transmission during the first 
channel use, A's transmission during the second channel use, 
B's transmission during the second channel use. We now 
represent these clusters by a Latin Square of side 16, with 
(xa 1 , xa 2 ) along the rows, and (xb 17 xb 2 ) along the columns. 
The {{xa 1 , xa 2 ), {xbx , x b 2 )) entry of the Latin Square as 
dictated by the clusters above are as follows: 



£i := {((0,0), (0,0)), ((0,1), (0,2 



C 2 



c 5 



c 6 



C 7 



((1,0), 
((2,0), 
((3,0), 

= {((0,0) 
((1,0), 
((2,0), 
((3,0), 

= {((0,0) 
((1,0), 
((2,0), 
((3,0), 

= {((0,0) 
((1,0), 
((2,0), 
((3,0), 

= {((0,0) 
((1,0), 
((2,0), 
((3,0), 

= {((0,0) 
((1,0), 
((2,0), 
((3,0), 

= {((0,0) 
((1,0), 
((2,0), 
((3,0), 



2,0)), ((1,1), (2, 2) 
1,0)), ((2, 1),(1,2) 

3.0) ), ((3, 1),(3,2) 

(0,3)), ((0,1), (0,1 

2, 3)), ((1,1), (2,1) 
1,3)), ((2, 1),(1,1) 

3, 3)), ((3, 1),(3,1) 

(0, 1)), ((0,1), (0,3 

2.1) ), ((1,1), (2, 3) 
1,1)), ((2, 1),(1,3) 

3.1) ), ((3, 1),(3,3) 

(0,2)), ((0,1), (0,0 

2, 2)), ((1,1), (2,0) 

1.2) ), ((2, 1),(1,0) 

3, 2)), ((3, 1),(3,0) 

(3.0) ), ((0,1), (3, 2 
1,0)), ((1,1), (1,2) 

2.0) ), ((2, 1),(2,2) 
0,0)), ((3, 1),(0,2) 

(3, 3)), ((0,1), (3,1 

1.3) ), ((1,1), (1,1) 
2, 3)), ((2, 1),(2,1) 
0,3)), ((3, 1),(0, 1) 

(3.1) ), ((0,1), (3, 3 

1.1) ), ((1,1), (1,3) 
2,1)), ((2, 1),(2,3) 
0,1)), ((3, 1),(0,3) 



), ((0,2), (0,1 
,((1,2), 
,((2,2), 
,((3,2), 



),((0,2) 
,((1,2), 
,((2,2), 
,((3,2), 

),((0,2) 
,((1,2), 
,((2,2), 
,((3,2), 

),((0,2) 
,((1,2), 
,((2,2), 
,((3,2), 

),((0,2) 
,((1,2), 
,((2,2), 
,((3,2), 

),((0,2) 
,((1,2), 
,((2,2), 
,((3,2), 

),((0,2) 
,((1,2), 
,((2,2), 
,((3,2), 



2,1) 
1,1) 
3,1) 

(0,2 
2,2) 
1,2) 
3,2) 

(0,0 
2,0) 
1,0) 
3,0) 

(0,3 
2,3) 
1,3) 
3,3) 

(3,1 
1,1) 
2,1) 
0,1) 

(3,2 
1,2) 
2,2) 
0,2) 

(3,0 
1,0) 
2,0) 
0,0) 



(0,3), (0,3)), 
1,3), (2, 3)), 

2, 3), (1,3)), 

3, 3), (3, 3))} 

(0,3), (0,0)), 
1,3), (2,0)), 

2, 3), (1,0)), 

3, 3), (3,0))} 

(0,3), (0,2)), 
1,3), (2, 2)), 

2, 3), (1,2)), 

3, 3), (3, 2))} 

(0,3),(0, 1)), 
1,3), (2,1)), 

2, 3), (1,1)), 

3, 3), (3,1))} 

(0,3), (3, 3)), 
1,3), (1,3)), 

2, 3), (2, 3)), 

3, 3), (0,3))} 

(0,3), (3,0)), 
1,3), (1,0)), 

2, 3), (2,0)), 

3, 3), (0,0))} 

(0,3), (3, 2)), 
1,3), (1,2)), 

2, 3), (2, 2)), 

3, 3), (0,2))} 
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(0,0) 


(0,1) 


(0,2) 


(0,3) 


(1,0) 


(1,1) 


(1,2) 


(1,3) 


(2,0) 


(2,1) 


(2,2) 


(2,3) 


(3,0) 


(3,1) 


(3,2) 


(3,3) 


(0,0) 


£i 


£3 


£4 


£2 


£9 


£11 
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£13 
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£5 


£7 
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Fig. 3. Latin Square L representing the clustering at the relay for the case je j0 = j, obtained using Cartesian Cluster 
Product, with the 4-PSK symbols that A(B) sent in the first and second channel use along the rows(columns) 



£ 8 := {((0, 0), (3, 2)), ((0, 1), (3, 0)), ((0, 2), (3, 3)), ((0, 3), (3, 1)), 
((1, 0), (1, 2)), ((1, 1), (1, 0)), ((1, 2), (1, 3)), ((1, 3), (1, 1)), 
((2, 0), (2, 2)), ((2, 1), (2, 0)), ((2, 2), (2, 3)), ((2, 3), (2, 1)), 
((3, 0), (0, 2)), ((3, 1), (0, 0)), ((3, 2), (0, 3)), ((3, 3), (0, 1))} 



£9 := {((0, 0), (1, 0)), ((0, 1), (1, 2)), ((0, 2), (1, 1)), ((0, 3), (1, 3)), 
((1, 0), (3, 0)), ((1, 1), (3, 2)), ((1, 2), (3, 1)), ((1, 3), (3, 3)), 
((2, 0), (0, 0)), ((2, 1), (0, 2)), ((2, 2), (0, 1)), ((2, 3), (0, 3)), 
((3, 0), (2, 0)), ((3, 1), (2, 2)), ((3, 2), (2, 1)), ((3, 3), (2, 3))} 



£ 15 := {((0, 0), (2, 1)), ((0, 1), (2, 3)), ((0, 2), (2, 0)), ((0, 3), (2, 2)), 
((1, 0), (0, 1)), ((1, 1), (0, 3)), ((1, 2), (0, 0)), ((1, 3), (0, 2)), 
((2, 0), (3, 1)), ((2, 1), (3, 3)), ((2, 2), (3, 0)), ((2, 3), (3, 2)), 
((3, 0), (1, 1)), ((3, 1), (1, 3)), ((3, 2), (1, 0)), ((3, 3), (1, 2))} 

£16 := {((0, 0), (2, 2)), ((0, 1), (2, 0)), ((0, 2), (2, 3)), ((0, 3), (2, 1)), 
((1, 0), (0, 2)), ((1, 1), (0, 0)), ((1, 2), (0, 3)), ((1, 3), (0, 1)), 
((2, 0), (3, 2)), ((2, 1), (3, 0)), ((2, 2), (3, 3)), ((2, 3), (3, 1)), 
((3, 0), (1, 2)), ((3, 1), (1, 0)), ((3, 2), (1, 3)), ((3, 3), (1, 1))} 



£ 10 := {((0, 0), (1, 3)), ((0, 1), (1, 1)), ((0, 2), (1, 2)), ((0, 3), (1, 0)), 
((1, 0), (3, 3)), ((1, 1), (3, 1)), ((1, 2), (3, 2)), ((1, 3), (3, 0)), 
((2, 0), (0, 3)), ((2, 1), (0, 1)), ((2, 2), (0, 2)), ((2, 3), (0, 0)), 
((3, 0), (2, 3)), ((3, 1), (2, 1)), ((3, 2), (2, 2)), ((3, 3), (2, 0))} 

£n := {((0, 0), (1, 1)), ((0, 1), (1, 3)), ((0, 2), (1, 0)), ((0, 3), (1, 2)), 
((1, 0), (3, 1)), ((1, 1), (3, 3)), ((1, 2), (3, 0)), ((1, 3), (3, 2)), 
((2, 0), (0, 1)), ((2, 1), (0, 3)), ((2, 2), (0, 0)), ((2, 3), (0, 2)), 
((3, 0), (2, 1)), ((3, 1), (2, 3)), ((3, 2), (2, 0)), ((3, 3), (2, 2))} 
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Fig. 4. Latin Square representing the clustering C^'l with 
the symbol sent by A(B) along the rows(columns). 



£12 := {((0, 0), (1, 2)), ((0, 1), (1, 0)), ((0, 2), (1, 3)), ((0, 3), (1, 1)), 
((1, 0), (3, 2)), ((1, 1), (3, 0)), ((1, 2), (3, 3)), ((1, 3), (3, 1)), 
((2, 0), (0, 2)), ((2, 1), (0, 0)), ((2, 2), (0, 3)), ((2, 3), (0, 1)), 
((3, 0), (2, 2)), ((3, 1), (2, 0)), ((3, 2), (2, 3)), ((3, 3), (2, 1))} 



£13 := {((0, 0), (2, 0)), ((0, 1), (2, 2)), ((0, 2), (2, 1)), ((0, 3), (2, 3)), 
((1, 0), (0, 0)), ((1, 1), (0, 2)), ((1, 2), (0, 1)), ((1, 3), (0, 3)), 
((2, 0), (3, 0)), ((2, 1), (3, 2)), ((2, 2), (3, 1)), ((2, 3), (3, 3)), 
((3, 0), (1, 0)), ((3, 1), (1, 2)), ((3, 2), (1, 1)), ((3, 3), (1, 3))} 



£14 := {((0, 0), (2, 3)), ((0, 1), (2, 1)), ((0, 2), (2, 2)), ((0, 3), (2, 0)), 
((1, 0), (0, 3)), ((1, 1), (0, 1)), ((1, 2), (0, 2)), ((1, 3), (0, 0)), 
((2, 0), (3, 3)), ((2, 1), (3, 1)), ((2, 2), (3, 2)), ((2, 3), (3, 0)), 
((3,0), (1,3)), ((3,1), (1,1)), ((3,2), (1,2)), ((3,3), (1,0))} 



The resulting Latin Square representing the clusters denoted 
by say L, is shown in Fig. 3. This 16 x 16 array L can be 
divided into 16 blocks of 4 x 4 arrays. Let Lb = \Li,j\ where 
each Zrjj is a 4 x 4 array for i,j — 1,2,3,4 as shown in 
Fig. 3. Each Li j is in one-to-one correspondence with the 
Latin Square obtained in [8| for removing the singular fade 
state jeJ 8 = j for the two-way 2-stage relaying scenario 
representing the clustering C^l given in Fig |4(a)| Also let, 
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= L^2 



This makes Lb of the form shown in Fig. |4(b)| 
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Fig. 5. Latin Square representing the clustering at the relay for the case -feJ e = 0.5 + 0.5j, obtained using Cartesian 
Cluster Product, with the 4-PSK symbols that A(B) sent in the first and second channel use along the rows(columns) 



This makes the block matrix Lb also consistent with the 
Latin Square given in Fig. 4. The reason behind Lb and 
Lij being consistent with this Latin Square is as follows: 
each L, j corresponds to some fixed values of the symbols 
A and B send during the first channel use, with the symbols 
sent by A and B during second channel use varying along 
the rows and columns respectively. The Latin Square in 
Fig. 3 has been obtained by taking the Cartesian Product 
of the clustering for removing the fade state ^e^ 9 = j with 
itself. The Cartesian Product utilizes the clustering that is 
represented by the Latin Square given in Fig. 4, given by 
C^'l, for the case for both the first and second channel use 
in the MA phase, which makes each Li j i,j — 1, 2, 3, 4 and 
Lb in one-to-one correspondence with this Latin Square that 
represents the clustering C^. 

Case 2: je je lies on the circle of radius 1 / \f2. 

In this case, the Cartesian Product of cl 7e3 1 with itself 

r jo] 

consists of 25 clusters since the clustering C' 7e ' has 5 
clusters. We now give an example of this case. The remaining 
instances of this case can be obtained from this example as 
will be shown later in Section V, Lemma 6. 

Example 2: Consider the case when je^ 6 = 0.5 + 0.5j. The 
clustering c[°- 5 +°- 5 i] given in |8| that removes this fade state 
for the two-way 2-stage relaying scenario is given by: 

£[0.5+0.5,] ={l u l 2 J 3 , k J 5 }, 

where, (11) 

?x ={(0,1), (1,2), (2, 3)} 
l 2 = {(0,2), (1,3), (3,0)} 
Z 3 = {(0,3), (2,0), (3,1)} 
2 4 = {(1,0),(2,1),(3,2)} 
h = {(0,0), (1,1), (2, 2), (3, 3)}. 



The Cartesian Product of the above clustering given by 
p[o.5+o.5j] = { C {u,i 3 } j 1,^ = 1,2,3,4,5} contains exactly 
25 clusters. The clusters and the corresponding constraints 
for the Latin Square representing the clustering have been 
listed in the Appendix A. The Cartesian Product of the 
clustering £[0-5+o.5i] with itself, denoted by £>[0-5+0.5i] can 
be represented by the Latin Square given in Fig. 5. 
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Fig. 6. Latin Square I representing the clustering 
c [o.5+o.5j] witn tne symbol sen t by A(B) along the 

rows(columns) 

As explained in the previous example, let the 16 x 16 
Latin Square obtained as shown in Fig. 5 be denoted by 
say L' and L' B = [L' id ] with i,j = 1,2,3,4. Then both 
L' B and L\ j must be consistent with the Latin Square 
of side 4 given in Fig. 6, denoted by which represents 
the clustering c[°- 5 + - 5 i]. As can be seen in Fig. 5, I is 
repeated in each block L\ j with a possibly different set 
of five symbols amongst {£1, £ 2 , £25} denoting the 
five symbols {k, k, ^5} in each L- ^. More precisely, 
the blocks L\ 1 = L' 22 = £33 = L\ 4 are the same 
as /, with the symbols £21, £22, £25 replacing the 
symbols k,k,---,k respectively. Similarly, the blocks 
L\ 2 = L' 2 3 = L' 3 4 are the same as / with symbols 
Cx,L-2,... ,£5 replacing the symbols k,k,---,k respectively, 
the blocks L' x 3 = L' 2 4 = L'± 1 are the same as I with symbols 
£6, £7, £10 replacing the symbols k,k,---,k respectively, 
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Fig. 7. Latin Square representing the clustering at the relay for the case je je = 1 + j, obtained using Cartesian Cluster 
Product, with the 4-PSK symbols that A(B) sent in the first and second channel use along the rows(columns) 



the blocks L\ 4 = 
symbols Cxi, £12, ■■ 



V, 



3.1 
£15 



L' A n are the same 



I 



replacing the symbols Zi , Z2 , 



with 

,.,/ 5 



respectively, and the blocks L' 2 1 = £3 2 = £4 3 are the same 
as / with symbols £16, £17, £20 replacing the symbols 
{i,22, ...,{5 respectively. Thus, the array £' can be obtained 
using I by simply using a different set of five symbols to 
denote li,l2 > —,h for every set of blocks corresponding to a 
symbol amongst l\, I2, ^5 in I. We will illustrate this in the 
next example by obtaining the 16 x 16 Latin Square using 
the 4x4 Latin Square given in [8| for the case. 
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Fig. 8. Latin Square representing the clustering C [1+jl 
with the symbol sent by A(B) along the rows(columns) 



Case 3: ^e je lies on the circle of radius \/2. 

In this case, the Cartesian Product of C^ 1 ^ 1 with itself 
consists of 25 clusters since the clustering C^ 1 ^ 1 has 5 
clusters. An instance of this case is as follows. 

Example 3: Consider the case when 7e J0 = 1 + j. The 
clustering C[ 1+ -?] given in |8| that removes this fade state for 
the two-way 2-stage relaying scenario is given by: 

C [1+Jl ={h,h,h,k,h} 
where, 
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= {(0,1) 


(2,3) 


(3,0)} 
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= {(0,3) 


(1,0) 


(3,2)} 
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= {(1,2) 


(2,0) 


(3,1)} 


k 


= {(0,2) 


(1,3) 


(2,1)} 
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= {(0,0) 


(1,1) 


(2, 2), (3, 3)} 



This clustering can be represented by a Latin Square of 
side 4 denoted by V as shown in Fig. 8. 

The Cartesian Product of the above clustering given by 
= {CfeW j ij = 1,2,3,4,5} contains exactly 25 
clusters as given in the Appendix B. We represent these 



clusters by a Latin Square of side 16, with (xa ± , %a 2 ) 
along the rows, and {xb 1 ,xb 2 ) along the columns. The 
((xax , xa 2 ), ( x Bi , xb 2 )) entry of the Latin Square as dictated 
by the clusters above are also listed in the Appendix B. 

Let the 16 x 16 Latin Square that represents the clustering 
obtained as a Cartesian Product of C^+J] with itself be 
denoted by say £", and L" B = [L? •] with i,j = 1,2,3,4. 
Then each of L" B and L'l must be consistent with the Latin 
Square of side 4 given in Fig. 8 which represents the clustering 
C^+J], We denote l\,l2, ~-,k in £1.2 = £3,4 = £4.1 by 
£i,£2, —,£5, in Li t 4 = L 2 .i = £-4,3 by £5, £6, £10, in 
blocks L 2 ,3 = £3,1 = £3,2 by £i , £11, £15, in blocks 
£•1,3 = £-2,4 = £3,2 by £i 6 ,£i7, ...,£ 2 o and in blocks 
£1,1 = £2,2 = £3,3 = £-4,4 by £21, £22, ■•-,£25- Placing 
these blocks in accordance with I', the Cartesian Product of 
the clustering C^+J] with itself, denoted by XJ^+J] can be 
represented by the Latin Square given in Fig. 7. 

The Latin Square which removes the singular fade state 
-e~i e can be obtained by taking the transpose of the 
Latin Square which removes the singular fade state je^ . 
For example, the Latin Square which removes the singular 
fade state ^eJi can also be obtained by taking the 
transpose of the Latin Square which removes the singular 
fade state -4=e~ J 'i. The reason for this is as follows: 
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the case when the singular fade state is ie~ j6 can be (j,9'), can be obtained from the Latin Square L of order 

2 2A which removes the singular fade state (7, 8), where 
1 — fe|x, as follows: Cyclic shift the columns of 
each one the 2 2A Latin Squares L,.j,0 < i,j ' < 2 A — 1, k 
times to the left to get the Latin Square L'. Cyclic shift the 
columns of the Square L' g associated with L', k times to the 
left, to get the Square L" B associated with the Latin Square L" . 



equivalently viewed as the case when the singular fade states 
is je^ with the users A and B interchanged. Interchanging 9' 
the users is equivalent to taking transpose of the Latin Square. 



IV. Clusterings from Latin Square of Lower Size 

In this section, we deal with any 2 A -PSK constellation. 
In 0, it was shown that the singular fade states for the 
two-way 2-stage relaying scenario lie on circles centered at 
the origin. From Lemma 1, since the singular fade states 
for the ACF two-way relaying scenario are the same as that 
of the two-way 2-stage relaying scenario, it follows that the 
singular fade states for the ACF two-way relaying scenario 
lie on circles centered at the origin as well. In this section, it 
is shown that for each circle, it is enough if we obtain one 
Latin Square which removes a singular fade state on that 
circle. The Latin Squares which remove the other singular 
fade states on that circle can be obtained by some elementary 
operations on the Latin Square, which are described in the 
sequel. 

For a Latin Square L of order 2 2A , let Lij, <i,j < 2 A -1, 
denote the Latin Sub-square of order 2 A obtained by taking 
only the rows 2 A i to 2 A (i + 1) — 1 and only the columns 
2 x j to 2 A (j + 1) - 1 of L. For example, in Fig. [9] the Latin 
Sub-squares Lij,i,j g {0,1}, of order 2 corresponding to 
the Latin Square L of order 4 are shown. Let Lb denote 
the Square of order 2 A , associated with the Latin Square L 
of order 2 2A , with £;.j,0 < i,j < 2 A — 1, as its entries. 
For example the square Lb of order 2 associated to a Latin 
Square L of order 4 is as shown in Fig. [10] 



-^0,0 
1 



1 



Fig. 9. Obtaining Latin Squares Lij's from the Latin Square L. 



£0,0 






Li,i 



Fig. 1 0. The Square L B order 2 corresponding to a Latin 
Square L of order 4. 

Lemma 3: Consider the two-way ACF relaying with 
2 A -PSK signal set used at nodes A and B. The Latin Square 
L" of order 2 2A which removes the singular fade state 



Proof: For the singular 

{((xa 1 ,xa 2 ),(x Bi ,x B2 )) , (jx' Al 
singularity removal constraint 



fade state 



i.e 



Br 



Ml 



7e 



be 



let 

a 



XA 2 



IB 2 -I 



x Bi x B t 

From (fL2l . it follows that 



_ jk2n _ jk2n 

(x Al ,x A2 ),(x Bl e 2* ,x B2 e 2* ) 



(12) 



B 2 



is a singularity removal constraint for the singular fade state 
(7,0'), where 9' — 9 = fc|f. In other words, the rotation 
in the je 1 * 8 plane by an angle 9 can be viewed equivalently 
as a rotation of the constellations used by B during the MA 
phases by an angle 9' — 8. Note that the columns of the Latin 
Square L which removes the singular fade states 7e 3 are 
indexed by the symbols (xb 1 ,xb 2 ) transmitted by B during 
the two MA phases. Rotating the signal set used by B during 
the second MA phase by an angle j s equivalent to shifting 
the columns of the Latin Sub-squares Lij k times to the left. 
Similarly, rotating the signal set used by B during the first 
MA phase by an angle ^ I is equivalent to cyclic shifting the 
columns of the square Lb, k times to the left. This completes 
the proof. 



For example consider the Latin Square L in Fig. |1 l(a)| 
which removes the singular fade state (7 = 1,9 = 0). The 
Latin Square L" which removes the singular fade state 
(7 = 1,9 = can be obtained from L as follows: The 
columns of the Latin Sub-squares Li j,0 < i,j < 3 are 
cyclically shifted once to the left, to obtain the Latin Square 
L' shown in Fig. |1 l(b)| The columns of the Square L' B 
associated with the Latin Square L' are cyclically shifted once 
to the left, to obtain the Square L" B . The Latin Square L" 
associated with the Square L" B which removes the singular 
fade state (7 = 1,9 = |) is shown in Fig. |ll(c)| 



For the ACF two-way relaying scenario with 4-PSK signal 
set used at the end nodes, the twelve singular fade states 
he on three circles with radii 1, -^=, and The Latin 
Squares which remove all the twelve singular fade states can 
be obtained from the three Latin Squares which remove the 
singular fade states j, 0.5 + 0.5j and 1 + j given in Fig. 3, 
Fig. 5 and Fig. 7 respectively. 
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V. Direct Clustering 



(a) The Latin Square L that removes the singular fade state (7 : 
1,0 = 0) 
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(b) The Latin Square V obtained form the Latin Square L using 
the procedure described in Lemma 4 
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(c) The Latin Square L" that removes the singular fade state 
(7 = l,0=f) 

Fig. 11. Construction of the Latin Square L" which 
removes (7 = 1,8 = tt/2) from the Latin Square L which 
removes (7 = 1,6 = 0) 



Algorithm 1: Obtaining the 16 x 16 Latin Square from 
the 16 x 16 array constrained using Singularity Removal 
Constraints 

Input: The constrained 16 x 16 array 

Output: A Latin Square representing the clustering map 
at the relay 

1 Start with the constrained 16 x 16 array 

2 Initialize all empty cells of the array to 

3 for 1 < i < 16 do 
4 
5 

6 

7 



9 

10 
11 
12 
13 

14 end 



for 1 < j < 16 do 

if cell of the array is NULL then 
Initialize c=l 

if C c does not occur in the i th row or the j 
column of the array then 

replace at cell of the array with 

else 

j c=c+l; 
end 

end 

end 



l.h 



Recall that there are three classes of singular fade states 
depending on the radius of the circle it lies on (Case 1:) je^ 
lies on the unit circle, (Case 2:) je^ e lies on the circle of radius 
l/i/2 and (Case 3:) -fe^ 6 lies on the circle of radius The 
number of clusters in the clustering utilized by relay node 
R during BC phase obtained using Cartesian Product in the 
three cases is 16, 25 and 25 respectively. It is observed that, 
if instead of taking the Cartesian Product of the clusterings 
given in (8j, the Cartesian Product of the Singularity Removal 
Constraints corresponding to each fade state are used to 
fill a 16 x 16 array, and the resulting incomplete array so 
obtained is completed using Algorithm [1] so as to form a 
Latin Square of side 16, then the number of clusters of the 
resulting clustering corresponding the to this Latin Square can 
be reduced from 25 to a lesser number in both Case 2 and 
Case 3. We call this the Direct Clustering. We now explain 
this clustering with the help of examples in the second and 
the third case only, since for the first case, the minimum 
number of clusters required, i.e., 16, can be achieved using 
Cartesian Product Clustering as shown in Section III. 



Case 2: je jd lies on the circle of radius 1/^/2. 

In this case, there are a total of 80 singularity removal 
constraints as shown in the following lemma. 

Lemma 4: When je^ 6 lies on the circle of radius 1/V2, 
there are a total of 80 singularity removal constraints. 

Proof: Let the singularity removal constraints for the two- 
way ACF relaying be the set {Ci, C2, C s }. Let, for t = 
l,2,...,s, 

Ct = {({xi k ,yi k ),(x2 k ,y2 k )) I xi k ,yi k ,X2 k ,y2 k G 5, k = 1,2, ...n t } . 
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Then, for 1 < fcijfcj < n t, 

x i kl + ie J0 y lki = x lk2 + je 30 y lk2 and 

x 2kl + ie >6 V2 kl = x 2k2 + ie je y 2k2 ■ 

Since in the case of two-way ACF relaying, the user nodes 
A and B transmit twice to the relay node R, these constraints 
for the ACF relaying can be obtained by taking the Cartesian 
Product of all sets of the form 

{(x Al ,x Bl ) | xa h +je je xB, 1 =x Al2 +je j9 x Bl2 VZi,Z 2 

where x All , x Bll , x Al2 , x Bl2 G S} . 
These sets can be of two types: 

1) The singularity removal constraints corresponding to the 
fade state 7e- 7 for the two-way 2-stage relaying as given 
in 0. Let us denote the 2-stage singularity removal 
constraints as: 

h = {(zii.S/uMa^ViJ} for i = 1,2,3,4. 

2) The sets of the form {(x A ,x B ), (x A ,x B )} for 
(x A ,x B ) € S 2 where (x A ,x B ) £ k V i = 1,2,3,4, 
since x A + ^e^xs = x A + je^xs- The (x a ,xb) for 
which (x A ,x B ) € U for some i = 1,2,3,4,5 are not 
considered in this category, as it already occurs in some 
set of the above category. 

The Cartesian Products of these sets amongst themselves 
must be the singularity removal constraints for the ACF 
relaying. Now, the constraint sets so obtained are also of two 
types: 

1) For i,j€ {1,2,3,4}, 

h x lj ={({x ll ,yi 1 ),{x' n ,y' n )),({x' il ,y' il ),(x jl ,y jl )), 
(On ,»<!), (x jl )), > Vix)) { x 'n . y'h ))} . 

These singularity removal constraints account for 16 
of the total number of constraints. 

2) For i € {1,2,3,4} and the four m k '■= 
((x Ak ,x Bk ),(xA k ,x Bk )) for k € {1,2, ...,8} that sat- 
isfy (x Ak ,x Bk ) i l 3 Vj G {1,2,3,4}; 

kxm k = {{{x ix ,y ix ), (xA h ,x Bk )), ((ai^.y-J, (x Ak ,x Bk ))} 
m k xli= {{(x Ak ,x Bk ), (xi^yij), ((x Ak ,x Bk ), (a^.y^))} 

These singularity removal constraints account for 
remaining 64 constraints. 

Thus, the set of singularity removal constraints for two-way 
ACF relaying becomes, 

{k xlj | i, j = 1, 2, 3, 4} U {li X m k | i = 1,2,3,4, k = 1, 2, 8} 
U {m fc X Z, | i = 1,2,3,4, k = 1, 2, 8} , 

where the subset {k x lj | i,j = 1,2,3,4} 

contains 16 constraints, and the subsets 

{li x m k \i = 1, 2, 3, 4, fc = 1, 2, 8} and 

{m k x Zj j t = 1,2,3,4, k = 1,2, ...,8} contain 32 

constraints each, which amount to a total of 80 singularity 



removal constraints. 



The 16 x 16 Latin Square representing these constraints can 
be completed using 20 symbols, as we show in the following 
example. 

Example 4: Consider the case for which je J = 
—0.5 + 0.5j. The singularity removal constraints for 
the case je^ 9 — —0.5 + 0.5j in two-way 2-stage relaying as 
given in [8 1 are: 

{(0,0), (1,3)}, {(1,1), (3, 2)}, {(0,1), (2, 2)} and {(2, 0), (3, 3)} . 

As a result, the singularity removal constraints for the 
two-way ACF relaying are as shown in Fig. 12 and Fig. 13. 
The clusters as shown in the third column of the table, are 
chosen such that each cluster satisfies the mutually exclusive 
laws given by (|5j and ©. 

The constraints can be represented using 17 symbols. In 
order to complete the latin square, we use Algorithm 1. A 
total of 20 symbols suffice in completing the array. Fig. 14 
represents the Latin Square representing the clustering at 
the relay for the case ^e> B = -0.5 + 0.5j, with the 4-PSK 
symbols that A(B) sent in the first and second channel use 
along the rows(columns). 

Case 3: je jd lies on the circle of radius 

In this case, there are a total of 80 singularity removal 
constraints as stated in the following lemma: 

Lemma 5: When ^e? s lies on the circle of radius V2, there 
are a total of 80 singularity removal constraints. 

We omit the proof of this lemma, as it is the same as that of 
Lemma 3. 

The resulting constrained 16 x 16 array can be uniquely 
completed with 20 symbols as shown in the following 
example: 

Example 5: Consider the case for which je^ e = — 1 + j. 
The singularity removal constraints for the case ^eJ = 
in two-way 2-stage relaying as given in 10 are: 

{(0,0), (3, 2)}, {(0,1), (3, 3)}, {(1,1), (2,0)} and {(1, 3), (2, 2)} . 

As a result, the constraints representing the singularity 
removal constraints for the two-way ACF relaying can be 
represented by the bold letters in the 16 x 16 Latin Square. 
The constraints can be represented using 18 symbols. In 
order to complete the latin square using Algorithm 1, a total 
of 20 symbols suffice. Fig. 15 represents the clustering. The 
constraints representing the singularity removal constraints 
for this example are given in Appendix C. 





Singularity Removal Constraints for 7e jy — —0.5 + 0.5j 


Latin Square Constraints for ■ye jU — —0.5 + 0.5j 
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1)), ((2, 3), (2, 2)), ((2, 3), (2, 2)) 


((0, 1) 


(1,1))} 


C 6 


(10) 


{((0, 1), (2, 0)) 


((2, 2), (3, 3)), ((2, 2), (3, 3)), ((0, 1), (2, 0))} 
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£11 


(28) 


{((1, 2) 


(0, 0)) 


((1, 2) 


(1,3))} 


{ ((1 , 0) , (2 
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{ ((1 , 0) , (2 


1)) 


((1, 2), (2, 2))} 




(45) 


{((2, 1), (0, 1)), ((2, 1), (2,2))} 


{ ((2, 0) , (1 


1)) 


((2, 2), (1, 2))} 


£14 


(46) 


{((2, 3), (0, 1)), ((2, 3), (2, 2))} 




1)) 


((2, 2), (3, 2))} 




(47) 


{((3,0), (0,1)), ((3,0), (2, 2))} 


{((3, 0) , (0 


1)) 


((3, 2), (0, 2))} 


£1 


(48) 


{((3, 1), (0, 1)), ((3, 1), (2, 2))} 


{((3, 0) , (1 


1)) 


((3, 2), (1, 2))} 


£11 


(49) 


{((1, 1), (0, 2)), ((3, 2), (0, 2))} 


{ ((1 , 0) , (1 


2)) 


((3,0), (2, 2))} 


£4 


(50) 


{((1, 1), (0,3)), ((3,2), (0,3))} 


{ ((1 , 0) , (1 


3)) 


((3,0), (2,3))} 




(51) 


{((1, 1), (1, 0)), ((3, 2), (1, 0))} 


{((1, 1), (1 


0)) 


((3, 1), (2,0))} 


£1.3 


(52) 


{((1, 1), (1, 2)), ((3, 2), (1, 2))} 


{((1, 1), (1 


2)) 


((3, 1), (2, 2))} 


£9 


(53) 


{((1, 1), (2, 1)), ((3, 2), (2, 1))} 


{ ((1 , 2) , (1 


1)) 


((3, 2), (2, 1))} 


£15 


(54) 


{((1, 1), (2, 3)), ((3, 2), (2, 3))} 


{ ((1 , 2) , (1 


3)) 


((3, 2), (2, 3))} 


£2 


(55) 


{((1,1), (3,0)), ((3, 2), (3,0))} 


{((1, 3) , (1 


0)) 


((3, 3), (2, 0))} 




(56) 


{((1, 1), (3, 1)), ((3,2), (3, 1))} 


{ ((1 , 3) , (1 


1)) 


((3, 3), (2, 1))} 


£1 


(57) 


{((0, 2), (1, 1)), ((0, 2), (3, 2))} 


{ ((0, 1) , (2 


1)),((0,3(,(2,2))} 




(58) 


{((0,3), (1, 1)), ((0,3), (3,2))} 


i 1 1 f\ / q 
H(U, i ) , yo 


1)), ((0,3), (3,2))} 


£12 


(59) 


{((1,0), (1, 1)), ((1,0), (3,2))} 


j 1 1 ~\ i in 
\yyi-, ±), y\) 


1)), ((1,3), (0,2))} 




(60) 


{((1,2), (1, 1)), ((1,2), (3,2))} 


J / 1 1 ~\\ lo 


1)), ((1,3), (2,2))} 


£17 


(61) 


{((2, 1), (1, 1)), ((2, 1), (3, 2))} 


J 1 1 O ~\\ 1 ~\ 

\\\*"> \*- 


1)), ((2, 3), (1,2))} 




(62) 


{((2,3) 


(1,1)) 


((2,3), (3,2))} 


{ ((2, 1) , (3 


1)) 


((2, 3), (3, 2))} 


£5 


(63) 


{((3,0) 


(1,1)) 


((3,0) 


(3,2))} 


{((3, 1) , (0 


1)) 


((3, 3), (0,2))} 


£n 


(64) 


{((3, 1) 


(1,1)) 


((3, 1) 


(3, 2))} 


{ ((3, 1) , (1 


1)) 


((3. 3), (1,2))} 


£12 


(65) 


{((2, 0) 


(0, 2)) 


((3, 3) 


(0, 2))} 


{ ((2, 0) , (0 


2)) 


((3.0), (3. 2))} 




(66) 


{((2, 0) 


(0,3)) 


((3,3) 


(0,3))} 


{ ((2, 0) , (0 


3)) 


((3,0), (3,3))} 


£e 


(67) 


{((2, 0) 


(1,0)) 


((3,3) 


(1,0))} 


{ ((2, 1) , (0 


0)) 


((3, 1), (3, 0))} 




(68) 


{((2, 0) 


(1,2)) 


((3,3) 


(1,2))} 


{ ((2, 1) , (0 


2)) 


((3, 1), (3, 2))} 


£15 


(69) 


{((2, 0) 


(2,1)) 


((3,3) 


(2,1))} 


{ ((2, 2) , (0 


1)) 


((3. 2), (3. 1))} 




(70) 


{((2, 0) 


(2, 3)) 


((3, 3) 


(2, 3))} 


{ ((2, 2) , (3 


3)) 


((3, 2), (3, 3))} 


£.3 


(71) 


{((2, 0) 


(3, 0)) 


((3, 3) 


(3,0))} 


{((2, 3) , (0 


0)) 


((3,3), (3,0))} 


£1.3 


(72) 


{((2,0), (3, 1)), ((3,3), (3, 1))} 


{ ((2, 3) , (0 


1)), ((3,3), (3, 1))} 




(73) 


{((0,2), (2,0)), ((0,2), (3,3))} 


-! ( ( 2) (2 


0)), ((0,3), (2,3))} 


£17 


(74) 


{((0,3), (2,0)), ((0,3), (3, 3))} 


{((0, 2), (3 


0)), ((0,3), (3,3))} 




(75) 


{((1,0), (2,0)), ((1,0), (3, 3))} 


{((1, 2), (0 


0)), ((1, 3), (0, :))} 


£12 


(76) 


{((1, 2), (2, 0)), ((1, 2), (3, 3))} 


{((1, 2), (2 


0)), ((1, 3), (2, 3))} 


£l6 


(77) 


{((2, 1), (2, 0)), ((2, 1), (3, 3))} 


{((2, 2), (1 


0)), ((2,3), (1,3))} 


£12 


(78) 


{((2, 3), (2,0)), ((2, 3), (3, 3))} 


{((2, 2), (3 


0)), ((2,3), (3,3))} 


£17 


(79) 


{((3,0), (2,0)), ((3,0), (3,3))} 


{((3, 2), (0 


0)), ((3,3), (0,3))} 


£17 


(80) 


{((3, 1), (2,0)), ((3, 1), (3,3))} 


{((3, 2), (1 


0)), ((3,3), (1,3))} 


£ 8 



Fig. 12. Singularity Removal Constraints Constraints for fe^ 9 = —0.5 + 0.5j 
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Fig. 13. Latin Square representing the clustering at the relay for the case ^e 30 obtained using Direct Clustering, with 
the 4-PSK symbols that A(B) sent in the first and second channel use along the rows(columns) 
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Fig. 14. Latin Square representing the clustering at the relay for the case jeJ e = -1 + j, with the 4-PSK symbols that 
A(B) sent in the first and second channel use along the rows(columns) 



VI. QUANTIZATION OF THE COMPLEX FADE 
STATE PLANE 

In practice, je^ 6 can take any value in the complex plane 
(it takes a value equal to one of the singular fade states 
with zero probability). As explained in Section II, one of 
the Latin Squares obtained, which remove the singular fade 
states needs to be chosen, depending on the value of je° . 
For a 7e je which is not a singular fade state, among all the 
Latin Squares which remove the singular fade states, the 
Latin Square C h ,h £ T~L which has the maximum value of 
the minimum cluster distance at je^ is chosen. In other 
words, for a given je^ ^ W, the clustering is chosen to 
be the one which removes the singular fade state h £ T-L 
which maximizes the metric d? min (C , 7e je ) given in (O. 



In this way, the je* — plane is quantized into \%\ point set, 
depending on which one of the obtained Latin Squares is 
chosen. 

For (x A ,x B ) j= { x 'a: x 'b) € £ 2 > let 
T>(j,8,xa,xb,x' a ,x' b ) be defined as, 

"D(l,0,x A ,XB,x' A ,x' B ) = \(x a - x' a ) + -)e ]0 {x B ~ x' B )\. 

(13) 

In the following lemma, it is shown that for a given 7e je , 
choosing the clustering C h £ C n , where h £ H, that maxi- 
mizes d? min (C h ,7e je ) given in (O, is the same as choosing 

the clustering CL XB x bJ ; where (xa,xb) 7^ ( x a> x 'b) e ^ 2 > 
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that minimizes the simpler metric given in ( fT3] l 
6: If the complex fa 

Cl x b- x bJ g c H are such that, 



Lemma 6: If the complex fade state je 3 and the clustering 

r *a-"a | 

L X B-*B. 



arg mm _ ^ T>(~y, 9, as a, %b , x' A , x' B ) = (x A , x B , x' A , x' B ), 



(xa,Xb)^(x' A} x' b )£S 2 



then 



G H, maximizes the metric d? min (C ,7e J ) 
given in (|9), among all h G "H. 

Proof: The squared minimum distance of the effective 
constellation at the relay d m i n (-fe 3 ) is given by ©. 

Let /i, /2 and /a be functions of 7 e :,e defined as in < TT~4t - 
( [ToT l given in the next page. We have, 

Ci„(7e iS ) = min{A( 7 e^), / 2 ( 7 e^), / 3 ( 7 e' 9 )}. 
From ( fT4b and ( TBI ), it follows that 



/i(7e j8 ) = /2(7e j9 )= min |( ¥ - x' B ) + ye> e (x B - ^ B )\ 

(xA,x B )i=(x A ,x B )!iS 2 



From ( [ToT l, it can be seen that, 



h{ie° 6 ) > min Kxa-^^^Vb-Zb)! 2 - 



Hence, we have 

d 2 ■ (~/e je ) 



min _ \(x A - x' B ) + ^e j0 {x B - x' B )\ 2 . 



{x A ,x B )j^(x' A ,x' B )£S 2 



Since, 



arg min 2? (7, 6>, ia, xb, x' a , x' b ) = (x A , x B , x' A , x' B ), 



we have, d 2 mm {ie 39 ) = \{x A - x' B ) + je je (x B - x 



/ \|2 



For the clustering C L XB " B J which removes the singular 
fade state — XA x f , the minimum cluster distance is greater 

XB Xg o 

than d m i n (je 38 ), while for all other clusterings in the set C u , 
it is equal to d m i n {^e 3 ). This completes the proof. 

■ 

The decision criterion in Lemma [6] based on which R 
chooses one of the Latin Squares obtained, is the same as 
the decision criterion for the two-way 2-stage relaying in 
Q- Hence, the quantization of the complex fade state plane 
for the ACF relaying is same as that of the two-way 2-stage 
relaying obtained in Q. 



VII. SIMULATION RESULTS 

The simulation results presented are for the case when 
Ha, Hb, H' a and H' B are distributed according to Rayleigh 
distribution, with the variances of all the fading links equal to 
dB. It is assumed that the AWGN noises at the three nodes 
are of variance dB. By SNR, we mean the average energies 
of the signal set used at the three nodes A, B and R, which 
are assumed to be equal. The frame length of a transmission 
is taken to be 256 bits. 



Consider the case when 4-PSK signal set is used at A and B. 
Fig. [TBI shows the SNR vs end-to-end sum throughput curves 
for the following schemes: Closest-Neighbour Clustering 
(CNC) Algorithm based scheme for the two-way 2-stage 
relaying proposed in 0, the Scheme based on Latin Squares 
for two-way 2-stage relaying proposed in (8), the scheme 
in which XOR network code is used irrespective of the 
channel condition, the Cartesian Product based scheme for 
ACF (ACF-CP) relaying and the Direct Clustering based 
scheme for ACF (ACF-DC) relaying. It can be seen from 
Fig. [TBJthat the schemes based on the ACF relaying perform 
better than the schemes based on 2-stage relaying at high 
SNR. From Fig. Q3] it follows that when the SNR is greater 
than 42 dB, the ACF-DC scheme outperforms all other 
schemes. The maximum throughput achieved by the ACF 
relaying schemes is 8/3 bits/s/Hz, whereas it is 2 bits/s/Hz 
for the 2-stage two-way relaying schemes. Also, as seen 
from Fig. Q3] the ACF-DC scheme performs better than the 
ACF-CP scheme. The reason for this is that the maximum 
cardinality of the signal set used during the BC phase is 25 for 
the ACF-CP scheme whereas it is 20 for the ACF-DC scheme. 

Consider the case when 8-PSK signal set is used at A and 
B. It was shown in [8| that for the two-way 2-stage relaying 
scheme with 8-PSK signal set, all the clusterings which 
the remove the singular fade states have exactly 8 clusters. 
Hence, the ACF-CP scheme, in which the clusterings are 
obtained by taking the Cartesian Product of the clusterings 
corresponding to the two-way 2-stage relaying scheme, have 
exactly 64 clusters (note that 64 is the minimum number of 
clusters required for conveying 6 information bits). Since 
the Cartesian Product itself results in the minimum number 
of clusters, the ACF-DC scheme is not considered for this 
case. Fig. [16] shows the SNR vs end-to-end sum throughput 
curves for the different scheme. Similar to the 4-PSK case, at 
high SNR, the ACF-CP scheme provides a larger throughput 
than the 2-stage relaying schemes. The maximum throughput 
achieved by the ACF-CP scheme is 4 bits/s/Hz, whereas it is 
3 bits/s/Hz for the 2-stage relaying schemes. 



VIII. Conclusion 

We proposed a scheme based on the ACF protocol for 
two-way relaying that utilizes totally three channel uses of 
the wireless two-way relaying channel unlike the 2-stage 
protocol that uses four channel uses, assuming that the users 
A and B transmit points from the same 4-PSK constellation. 
The network codes used at the relay during the Broadcast 
Phase were obtained using two methods: by taking the 
Cartesian Product of the clusterings proposed in [8] for 
two user 2-stage case and by completing the Latin Square 
filled partially with the singularity removal constraints for 
a given fade state. Using the second method called Direct 
Clustering, the maximum size of the resulting constellation 
used by the relay node R in the BC phase was reduced to 
20, as compared to the Cartesian Product based approach 
which results in the constellation size being 25 for these 
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flife 3 ") = min \\{xa 1 - x' A ) + -ye 3e {x Bl - x' )| 2 + |(x,4 - x' A ) + 7e :,( '(x B - x' )| 2 !• (14) 

C-Al ffl '^'A! fflj^Aj ^-B 2 ) = <^ 2 ^Bg ) * __> 

/2(7e J ") = , min I \{x A - x' A ) + -ye 3e {x Bl - x' )| 2 + \{x A - x' A ) + 7e 3,) (x B - x' B )| 2 !• (15) 
,- Bl )=(^! .* Bl >.(^ 2 ^b 2 )^(^ 2 ^'b 2 ) * * I 2 __> 

hiie 3 ") = min { |(x Al - x' A ) + 7e 3 "(x B - x' )| 2 + |(xa 2 - x' A ) + - 1 e 3 "(x B - x' B )| 2 !• (16) 

(-a 1 .-b 1 )^(-' Ai .-'b 1 ).(-a 2 .-b 2 )^<^ 2 .- B2 ) 1 2 J 
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Fig. 15. SNR vs throughput curves for different schemes for 4-PSK signal set 



cases. Having obtained all the Latin Squares, the complex 
plane was quantized depending on which one of the obtained 
Latin Squares maximizes the minimum cluster distance. This 
quantization was shown to be the same as that achieved in Q 
for the two-way 2-stage relaying scenario. Simulation results 
showed that the ACF protocol based schemes outperform the 
schemes proposed in [4] and |8|, at high SNR. 
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Appendix A 

Clustering that removes the singular fade state 0.5 + 0.5j." 
The Cartesian Product of the clustering c[°- 5 + - 5 j] with itself 
denoted by X>[o.5+o.5j] = | C {Wj} \ = 1,2,3,4,5} con- 
tains exactly 25 clusters. The clusters in X>[°- 5 +°- 5 j] have been 
listed below. 
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The entries in the above constraints are of the form 
{{xa 1 tXb 1 ) 1 (xa 2 ,Xb 2 ))- I n order to represent these clusters 
by a Latin Square of side 16, with {xa 1 , xa 2 ) along the rows, 
and (xb 1 ,xb 2 ) along the columns each entry must be made of 
the form {{xa 1 , xa 2 ), (xb 1 ,xb 2 ))- Thus the constraints on the 
Latin Square as dictated by the clusters above are as follows: 
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(1,2)), ((2, 3), 
(2, 2)), ((3, 3), 

• (0,3)), ((1,0) 

(2.1) ), ((2, 1), 

(3. 2) ), ((3, 2), 

• (0,0)), ((1,0) 
(2, 2)), ((2,1), 
(3, 3)), ((3, 3), 

• (0,1)), ((1,0) 

(2. 3) ), ((2, 2), 
(3,0)), ((3, 3), 

• (0,2)), ((1,1) 

(2.0) ), ((2, 2), 

(3.1) ), ((3, 3), 

• (0,2)), ((0,3) 

(1.2) ), ((1,3), 
(2, 2)), ((2, 3), 
(3, 2)), ((3, 3), 



,(0,3)), 

(1,3)), 

(2,3))} 

,(1,1)), 

(2.2) ), 

(3.3) )} 
,(1,2)), 
(2,3)), 
(3,0))} 
,(1,3)), 

(2.0) ), 

(3. 1) )} 
,(1,0)), 

(2.1) ), 

(3.2) )} 
,(0,3)), 

(1.3) ), 
(2,3)), 
(3,3))} 



Appendix B 

Clustering that removes the singular fade state 1 + j: 
The Cartesian Product of the clustering C^ 1+ ^ with itself 
denoted by V^ 1+ ^ contains exactly 25 clusters. The clusters 
in are as follows. 



c {hM = {((0 , i), (0, 1)), ((0, 1), (2, 3)), ((0, 1), (3, 0)), 
((2, 3), (0,1)), ((2, 3), (2, 3)), ((2, 3), (3,0)) 
((3,0), (0,1)), ((3,0), (2, 3)), ((3,0), (3,0)),} 

COi.b} = {((0, 1), (0, 3)), ((0, 1), (1, 0)), ((0, 1), (3, 2)), 
((2, 3), (0,3)), ((2, 3), (1,0)), ((2, 3), (3, 2)) 
((3,0), (0,3)), ((3,0), (1,0)), ((3,0), (3, 2)),} 

COi.W = {((0, 1), (1, 2)), ((0, 1), (2, 0)), ((0, 1), (3, 1)), 
((2, 3), (1,2)), ((2, 3), (2,0)), ((2, 3), (3,1)) 
((3,0), (1,2)), ((3,0), (2,0)), ((3,0), (3,1)),} 

C^M = {((0 , 1), (0, 2)), ((0, 1), (1, 3)), ((0, 1), (2, 1)), 
((2, 3), (0,2)), ((2, 3), (1,3)), ((2, 3), (2,1)) 
((3,0), (0,2)), ((3,0), (1,3)), ((3,0), (2,1)),} 

C Oi,W = {((0 , 1), (0, 0)), ((0, 1), (1, 1)), ((0, 1), (2, 2)), ((0, 1), (3, 3)), 
((2, 3), (0, 0)), ((2, 3), (1, 1)), ((2, 3), (2, 2)), ((2, 3), (3, 3)), 
((3, 0), (0, 0)), ((3, 0), (1, 1)), ((3, 0), (2, 2)), ((3, 0), (3, 3))} 

CV'M = {((0, 3), (0, 1)), ((0, 3), (2, 3)), ((0, 3), (3, 0)), 
((1,0), (0,1)), ((1,0), (2, 3)), ((1,0), (3,0)) 
((3, 2), (0, 1)), ((3, 2), (2, 3)), ((3, 2), (3, 0)), } 



C {h 



C {h 



C {h 



C {h 



C {h 



c {i 4 



c {i 4 



C {i 4 



c {! 4 



c {i 4 



C {h 



C° 5 



3} 



5} 



3} 



3} 



(0,3), 

1,0 

3,2 

(0,: 



(0,3 
0,3) 
0,3) 
(1,2 
1,2) 
1,2) 
(0,2 
0,2) 
0,2) 
(0,0 
0,0) 
0,0) 
(0,1 

0,1) 
0,1) 

(0,3 
0,3) 
0,3) 

(1,2 
1,2) 
1,2) 
(0,2 
0,2) 
0,2) 
(0,0 
0,0) 
0,0) 
(0,1 
0,1) 
0,1) 
(0,3 
0,3) 
0,3) 
(1,2 
1,2) 
1,2) 
(0,2 
0,2) 
0,2) 
(0,0 
0,0) 
0,0) 

(0,1 
2,3) 
3,0) 
(0,3 
1,0) 
3,2) 

(1,2 
2,0) 
3,1) 
(0,2 
1,3) 
2,1) 



((0,3) 

(1,0) 

(3,2) 

((0,3) 

(1,0) 

(3,2) 

((0,3) 

(1,0) 

(3,2) 

((0,3) 

(1,0) 

(3,2) 

((1,2) 

(2,0) 

(3,1) 

((1,2) 

(2,0) 

(3,1) 

((1,2) 

(2,0) 

(3,1) 

((1,2) 

(2,0) 

(3,1) 

((1,2) 

(2,0) 

(3,1) 

((0,2) 

(1,3) 

(2,1) 

((0,2) 

(1,3) 

(2,1) 

((0,2) 

(1,3) 

(2,1) 

((0,2) 

(1,3) 

(2,1) 

((0,2) 

(1,3) 

(2,1) 

((0,0) 

(1,1) 
(3,3) 

((0,0) 

(1,1) 
(3,3) 

((0,0) 

(1,1) 
(3,3) 

((0,0) 

(1,1) 
(3,3) 



(1,0 
1,0 

1,0 

(2,< 
2,0 
2,0 

(1,: 
1,3 
1,3 

(1,: 
1, 1 
1, 1 

(2,: 
2,3 
2,3 

(i,( 
1,0 
1,0 

(2,< 
2,0 
2,0 

(1,1 
1,3 
1,3 

(1,: 
1,1 
1, 1 

(2,1 
2,3 
2,3 

(1,< 
1,0 
1,0 

(2,< 
2,0 
2,0 

(1,: 
1,3 
1,3 

(1,: 
1, 1 
1,1 

(2,: 
3,0 
0,1 

(i,( 
3,2 
0,3 

(2,( 
3, 1 
1,2 
1,3 
2, 1 
0,2 



M(0,3), 
,((1,0) 
, ((3,2) 
),((0,3) 
,((1,0) 
, ((3,2) 
),((0,3) 
,((1,0) 
, ((3,2) 
),((0,3) 
,((1,0) 
, ((3,2) 

),((!, 2) 
, ((2,0) 
,((3,1) 
),((!, 2) 
, ((2,0) 
,((3,1) 

),((!, 2) 

, ((2,0) 

,((3,1) 

),((!, 2) 

, ((2,0) 

,((3,1) 

),((!, 2) 

, ((2,0) 

,((3,1) 

),((0,2) 

,((1,3) 

,((2,1) 

),((0,2) 

,((1,3) 

,((2,1) 

),((0,2) 

,((1,3) 

,((2,1) 

),((0,2) 

,((1,3) 

,((2,1) 

),((0,2) 

,((1,3) 

,((2,1) 

),((0,0) 

, ((2,2) 

, ((3,3) 

),((0,0) 

, ((2,2) 

,((3,3) 

),((0,0) 

, ((2,2) 

, ((3,3) 

, ((0,0) 

, ((2,2) 
,((3,3) 



(3,2; 

3,2); 
3,2); 

(3,i; 

3,1); 
3,1); 
(2,1; 
2,1); 

2.1) ; 
(2,2; 

2.2) ; 
2,2); 
(3,0; 
3,0); 

3.0) ; 
(3,2; 
3,2); 
3,2); 

(3,1; 

3.1) ; 
3,1); 
(2,1; 
2,1): 

2.1) : 
(2,2; 

2.2) ; 
2,2); 
(3,0; 
3,0); 

3.0) ; 
(3,2; 
3,2); 
3,2); 
(3,1; 

3.1) : 
3, i): 
(2,1: 
2,1): 

2.1) : 
(2,2; 

2.2) ; 

2.2) ; 
(3,0; 
0,1): 

2.3) ; 
(3,2; 
0,3); 
1,0); 
(3,1; 

1.2) ; 

2.0) ; 

2.1) : 

0,2); 

1.3) : 



), 

,} 

), 

,} 

), 

,} 

),((<), 3), (3, 3)), 
,((1,0), (3,3)), 
,((3, 2), (3, 3))} 

), 

,} 
), 

,} 

), 

,} 

), 

,} 

), ((1,2), (3,3)), 
,((2,0), (3, 3)), 
,((3,1), (3, 3))} 

), 

,} 

), 

,} 

), 

,} 

), 



• } 



((0,2), (3, 3)), 



(1,3 
(2,1 

((1,1) 

(2,2 

(3,3 

((1,1) 

(2,2 

(3,3 

((1,1) 

(2,2 

(3,3 

(1, 1 

(2,2 

(3,3 



(3,3)), 
(3,3))} 

,(0,1)), 

(2,3)), 

(3,0))} 

,(0,3)), 

(1,0)), 

(3.2) )} 

,(1,2)), 

(2.0) ), 

(3.1) )} 
(0,2)), 

(1.3) ), 
(2,1))} 
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c {h,i 5 } = {((o, o), (0, 0)), ((0, 0), (1, 1)), ((0, 0), (2, 2)), ((0, 0), (3, 3)), 
((1, 1), (0, 0)), ((1, 1), (1, 1)), ((1, 1), (2, 2)), ((1, 1), (3, 3)), 
((2, 2), (0, 0)), ((2, 2), (1, 1)), ((2, 2), (2, 2)), ((2, 2), (3, 3)), 
((3, 3), (0, 0)), ((3, 3), (1, 1)), ((3, 3), (2, 2)), ((3, 3), (3, 3))} 



The constraints for the Latin Square representing the above 
clustering are as follows: 



(0,0), (1,1 
2,0), (3,1) 
3,0), (0, 1) 
(0,0), (1,3 
(3,3) 
(0,3) 
(0,1), (1,2 
,(3,2) 
,(0,2) 



2,0 
3,0 



2.1 
3. 1 



2,0 
3,0 



(0,0), (1,2 
,(3,2) 
,(0,2) 
(0,0), (1,0 
,(3,0) 
,(0,0) 



2,0 
3,0 



(0,0), (3,1 



1,0 
3,0 



(0,0), (3, 3 



1,0 
3,0 



1,1 
3. 1 



2,0 
3,0 



2,0 
3,0 



2.1 
3. 1 



2,0 
3,0 



(0, 1) 
(2,1) 



(0,3) 
(2,3) 



(0,1), (3, 2 
,(0,2) 
,(2,2) 
(0,0), (3, 2 
,(0,2) 
,(2,2) 



1,0 
3,0 
(0,0), (3,0 



1,0 
3,0 



,(0,0) 
,(2,0) 
(1,0), (2,1 



, (0, 1) 
,(1,1) 



(1,0), (2, 3 



,(0,3) 
,(1,3) 



(1,1), (2, 2 
, (0, 2) 
,(1,2) 
(1,0), (2, 2 
,(0,2) 
,(1,2) 



(1,0), (2,0 
2,0), (0,0) 
3,0), (1,0) 



(0,2 
2,2) 
3,2) 
(0,1 
2,1) 
3,1) 
(0,2 
2,2) 
3,2) 
(0,1 
2,1) 
3,1) 
(0,1 
2,1) 
3,1) 
(0,2 
1,2) 
3,2) 
(0,1 
1,1) 
3,1) 
(0,2 
1,2) 
3,2) 
(0,1 
1,1) 
3,1) 
(0,1 
1,1) 
3,1) 
(1,2 
2,2) 
3,2) 

(1,1 
2,1) 
3,1) 
(1,2 
2,2) 
3,2) 
(1,1 
2,1) 
3,1) 
(1,1 
2,1) 
3,1) 



,(1,3 
(3,3) 
(0,3) 
,(1,0 
(3,0) 
(0,0) 
,(1,0 
(3,0) 
(0,0) 
,(1,3 
(3,3) 
(0,3) 
,(1,1 
(3, 1) 
(0, 1) 
,(3,3 
(0,3) 
(2,3) 
,(3,0 
(0,0) 
(2,0) 
,(3,0 
(0,0) 
(2,0) 
,(3,3 
(0,3) 
(2,3) 
,(3,1 
(0, 1) 
(2, 1) 
,(2,3 
(0,3) 
(1,3) 
,(2,0 
(0,0) 
(1,0) 
,(2,0 
(0,0) 
(1,0) 
,(2,3 
(0,3) 
(1,3) 
,(2,1 
(0, 1) 
(1,1) 



(0,3), (1,0 
2, 3), (3,0) 
3,3), (0,0) 
(0,3), (1,2 
2, 3), (3, 2) 
3,3), (0,2) 
(0,3), (1,1 

2, 3), (3,1) 

3, 3), (0,1) 
(0, 2), (1,1 

2, 2), (3,1) 

3, 2), (0,1) 
(0, 2), (1,2 

2, 2), (3, 2) 

3.2) , (0,2) 
(0,3), (3,0 

1.3) , (0,0) 
3,3), (2,0) 
(0,3), (3, 2 
1,3), (0,2) 
3,3), (2,2) 
(0,3), (3,1 
1,3), (0,1) 

3, 3), (2,1) 
(0, 2), (3,1 
1,2), (0,1) 
3, 2), (2,1) 
(0, 2), (3, 2 
12), (0,2)) 

3.2) , (2,2) 
(1,3), (2,0 
2, 3), (0,0) 

3.3) , (1,0) 
(1,3), (2, 2 
2, 3), (0,2) 
3,3), (1,2) 
(1,3), (2,1 

2, 3), (0,1) 

3, 3), (1,1) 
(1,2), (2,1 

2, 2), (0,1) 

3, 2), (1,1) 
(1, 2), (2,2 
2,2), (0,2) 
3, 2), (1,2) 



). 
,} 

), ((0,3), (1,3)), 
,((2, 3), (3, 3)), 
,((3, 3), (0,3))} 
), 

,} 



,} 

),((0,3),(3,3)), 
((1,3), (0,3)), 
,((3, 3), (2, 3))} 
), 

,} 



,} 

),((!, 3), (2, 3)), 
,((2, 3), (0,3)), 
,((3, 3), (1,3))} 



£l6 


= {((0,0) 


,(2,1)), 


((0,2), (2,3 


), ((0,3), (2,0 


), 






((1,0), 


(3, 1)), ( 


(1,2) 


(3,3) 


,((1,3) 


(3,0) 








((2,0), 


(1, 1)), ( 


(2,2) 


(1, 3) 


, ((2,3) 


(1,0) 


,} 




Cl7 


= {((0,0) 


,(2,3)), 


((0,1) 


,(2,0 


),((0,3) 


,(2,2 


) 5 






((1,0), 


(3,3)), ( 


(1, 1) 


(3,0) 


,((1,3) 


(3,2) 








((2,0), 


(1,3)), ( 


(2, 1) 


(1,0) 


, ((2,3) 


(1,2) 


,} 




£is 


= {((0,1) 


,(2,2)), 


((0,2) 


,(2,0 


M(0,3) 


,(2,1 


) 1 






((1,1), 


(3,2)),( 


(1,2) 


(3,0) 


,((1,3) 


(3, 1) 








((2, 1), 


(1,2)), ( 


(2,2) 


(1,0) 


, ((2,3) 


(1, 1) 


,} 




£l9 


= {((0,0) 


,(2,2)), 


((0,1) 


,(2,3 


),((o, 2; 


,(2,1 


) ! 






((1,0), 


(3,2)),( 


(1, 1) 


(3,3) 


,((1,2) 


(3, 1) 








((2,0), 


(1,2)), ( 


(2, 1) 


(1, 3) 


, ((2,2) 


(1, 1) 


,} 




£20 


= {((0,0) 


,(2,0)), 


((0,1) 


,(2,1 


),((o, 2; 


,(2,2 


),((0,3),(2,3)), 




((1,0), 


(3,0)), ( 


(1, 1) 


(3, 1) 


,((1,2) 


(3,2) 


,((1,3) 


(3,3)), 




((2,0), 


(1,0)), ( 


(2, 1) 


(1, 1) 


, ((2,2) 


(1,2) 


,((2,3) 


(1,3))} 


£21 


= {((0,0) 


,(0,1)), 


((0,2) 


,(0,3 


),((0,3) 


,(0,0 


),((1,0) 


,(1,1)), 




((1,2), 


(1,3)), ( 


(1,3) 


(1,0) 


, ((2,0) 


(2,1) 


, ((2, 2) 


(2,3)), 




((2,3), 


(2,0)), ( 


(3,0) 


(3, 1) 


, ((3,2) 


(3,3) 


,((3,3) 


(3,0))} 


£22 


= {((0,0) 


,(0,3)), 


((0,1) 


,(0,0 


),((0,3) 


,(0,2 


),((1,0) 


,(1,3)), 




((1,1), 


(1,0)), ( 


(1,3) 


(1,2) 


, ((2,0) 


(2,3) 


,((2, 1) 


(2,0)), 




((2,3), 


(2,2)), ( 


(3,0) 


(3,3) 


, ((3, 1) 


(3,0) 


,((3,3) 


(3,2))} 


£23 


= {((0,1) 


,(0,2)), 


((0,2) 


,(0,0 


),((0,3] 


,(0,1 


),((!,!) 


,(1,2)), 




((1,2), 


(1,0)), ( 


(1, 3) 


(1, 1) 


,((2,1) 


(2,2) 


,((2,2) 


(2,0)), 




((2,3), 


(2, 1)), ( 


(3, 1) 


(3,2) 


, ((3,2) 


(3,0) 


,((3,3) 


(3,1))} 


£24 


= {((0,0) 


,(0,2)), 


((0,1) 


,(0,3 


),((o, 2; 


,(0,1 


),((1,0) 


,(1,2)), 




((1,1), 


(1,3)),( 


(1,2) 


(1,1) 


,((2,0) 


(2,2) 


,((2, 1) 


(2,3)), 




((2,2), 


(2, 1)), ( 


(3,0) 


(3,2) 


,((3,1) 


(3,3) 


, ((3, 2) 


(3,1))} 


£25 


= {((0,0) 


,(0,0)), 


((0,1) 


,(0,1 


),((o, 2; 


,(0,2 


),((0,3) 


,(0,3)), 




((1,0), 


(1,0)), ( 


(1,1) 


(1,1) 


,((1,2) 


(1,2) 


,((1,3) 


(1,3)), 




((2,0), 


(2,0)),( 


(2,1) 


(2,1) 


,((2,2) 


(2,2) 


,((2,3) 


(2,3)), 




((3,0), 


(3,0)), ( 


(3,1) 


(3,1) 


,((3,2) 


(3,2) 


,((3,3) 


(3,3))} 



The Cartesian Product of the clustering Cl 1+ ^ with itself, 
denoted by 2?[ 1+ ^ can be represented by the Latin Square 
given in Fig. 8. 



Appendix C 

Singularity removal constraints for the singular fade state 
-l+j: ' 

The singularity removal constraints for the singular fade state 
— l+ j are given in Fig. 17 in the next page. 





Singularity Removal Constraints for 7e jy — — 1+j 


Latin Square Constraints for ~ye Jti — — 1+ j 


Cluster 


Ml 


{((0, 0), (3, 2)), ((3, 2), (0, 0)), ((0, 0), (0, 0)), ((3, 2), (3, 2))} 


ilia ^l fTl 9')') 111 O 1 ! 19 nil 1(f) HI fTl 0)) 111 1) 19 91U 


r , 


C91 


{((0, 0), (0, 1)), ((3, 2), (3, 3)), ((3, 2), (0, 1)), ((0, 0), (3, 3))} 


iim ni i^ ^ 19 n 1 ! 19 ^)) no id s^h 


r 


w) 


{((0, 0), (1, 1)), ((3, 2), (2, 0)), ((3, 2), (1, 1)), ((0, 0), (2, 0))} 


iim i 1 ) in i 1 ) 1 ) (a 9} 19 (a 19 ^\\ no 9) la niu 


*-3 


f41 


{((0, 0), (1, 3)), ((3, 2), (2, 2)), ((3, 2), (1, 3)), ((0, 0), (2, 2))} 


illf) 1 ^ fTl (Ci 9} 19 9\\ (Ci 1 ^ 19 lit) 9) Id 9))\ 




rsi 


{((0, 1), (0,0)), ((3, 3), (3, 2)), ((3, 3), (0,0)), ((0, 1), (3,2))} 






f6i 


{((0, 1) 


(1, 1)) 


((3,3) 


(2, 0)), ((3, 3), (1, 1)), ((0, 1), (2, 0))} 


{((0, 1), (1, 1)), ((3, 2), (3, 0)), ((3, 1), (3, 1)). ((0.2). (1.0))} 


*— i 


\ 1 > 


{((0, 1) 


(0, 1)) 


((3,3) 


(3,3)), ((3,3), (0, 1)), ((0, 1), (3,3))} 


-!f(0 01 M 111 (73 31 f3 311 ff3 01 (3 1 ) ) 170 31 C\ 31H 




(R\ 
\ a ) 


{((0, 1) 


(!.3)) 


((3,3) 


(2, 2)), ((3, 3), (1, 3)), ((0, 1), (2, 2))} 


11 1 "^11 (71 91 f'l 911 f'f'l 11 fl 111 (fn 91 M 91H 


T„ 
*-2 




{((!. !) 


(1, 1)) 


((2,0) 


(2, 0)), ((2, 0), (1, 1)), ((1, 1), (2, 0))} 


/(■(■l 11 f 1 111 ^9 91 f'O Oil (79 11 f'O ill i7l 91 ( 1 fl11\ 

1111) v-L; -'-//; vv^j^/i v ' / / ' \ v ) / ) V ' //'V\ ' /'\ ' //j 
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